Charge renormalization and other exact coupling corrections in the dipolar efFective 
interaction in an electrolyte near a dielectric wall 
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C ' The aim of the paper is to study the renormalizations of the charge and of the screening length 

' O ' that appear in the large-distance behavior of the effective pairwise interaction w„„/ between two 

O I charges e^ and e„' in a dilute electrolyte solution, both along a dielectric wall and in the bulk. The 

electrolyte is described by the primitive model in the framework of classical statistical mechanics 

^ , and the electrostatic response of the wall is characterized by its dielectric constant. In Ref.|l| a 



fT^ ' graphic reorganization of resummed Mayer diagrammatics has been devised in order to exhibit the 

\0 ' general structure of the 1/y leading tail oi w^^i{x,x' ,y) for two charges located at distances x and 

o, 

^^ , x' from the wall and separated by a distance y along the wall. When all species have the same closest 



approach distance b to the wall, the coefficient of the 1/y tail is the product Dc{x)D^i{x') of two 
effective dipoles. Here we use the same graphic reorganization in order to systematically investigate 



' O ' the exponential large-distance behavior of w„„' in the bulk. (We show that the reorganization also 

O ■ enables one to derive the basic screening rules in both cases.) Then, in a regime of high dilution 

and weak coupling, the exact analytical corrections to the leading tail of w^^i , both in the bulk or 

^\ ' along the wall, are calculated at first order in the coupling parameter e and in the limit where b 

becomes negligible with respect to the Debye screening length, {e is proportional to the so-called 
plasma parameter.) The structure of corrections to the terms of order e is exhibited, and the scaling 
regime for the validity of the Debye limit is specified. In the vicinity of the wall, we use the density 
profiles calculated in Ref.|2| up to order e and the method devised in Ref.|3| for the determination 
of the corresponding correction in the auxiliary screened potential, which also appears in the linear- 



response theory. The first coupling correction to the effective dipole Da(x) is a function (not a 
mere exponential decay) determined by the nonuniformity of the density profiles as well as by three- 
and four-body screened interactions in w^^i . Though the effective screening length (beyond the 
Debye value) in the direction perpendicular to the wall is the same as in the bulk, the bare solvated 
charges are not renormalized by the same quantity as in the bulk, because of combined steric and 
electrostatic effects induced by the wall. 

PACS numbers: 05.20. Jj, 05.70.Np. 

Corresponding author : cornu@th.u-psud.fr. 

I. INTRODUCTION 



A. Issue at stake 

The paper is devoted to the large-distance behavior of the pairwise eflective interaction betvi^een two charges in an 
electrolyte solution, which is confined to the region a; > by a plane impenetrable dielectric wall. The electrolyte 
solution is described by the usual primitive model [j| with rig species of charges which interact via the Coulomb 
interaction. Every charged particle of species a is represented as a hard sphere - with diameter a^ - where the net 
bare solvated charge e„ = Z„e is concentrated at the center of the sphere, (e is the abolute value of the electron 
charge and Z„ may be positive or negative.) The solvent is handled with as a continuous medium of uniform dielectric 
constant esoiv The wall matter is characterized by a dielectric constant e„, 7^ egoiv, and the latter difference results 
into an electrostatic response of the wall to the moving charges in the electrolyte. Moreover, the excluded-volume 
sphere of every particle is assumed to be made of a material with the same dielectric constant as that of the solvent. 
(Therefore e = Csoiv when x > and e = e^ when a; < 0). In the framework of statistical mechanics, the effective 
pairwise interaction u;„„/(r,r') between two charges e„ and e„' located at positions r and r', respectively, is defined 

ri 

from the pair correlation function h^^r by (see e.g. Ref.[a|) 

1 -I- /i„„/ = exp(-/3K;„„/), (1) 
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where /? — l/k^T is the inverse temperature, in which k^ is Boltzmann constant and T is the absolute temperature. 
(w„q' is also called potential of mean force, while h^^> is known as the Ursell function.) In the vicinity of the wall, 
symmetries enforce that w„„/(r, r') — w^^i {x, x' ,y), where x and x' are the distances of r and r' from the wall and y 
is the norm of the projection y of r — r' onto the wall plane. Along the wall, contrary to the bulk case, «;„„' (x, x', y) 
does not decay exponentially fast: its leading behavior at large distances y takes a dipolar form f^^i{x,x')/y'^ , as 
a result of the deformation of screening clouds enforced by the presence of the wall (see Ref|^] for a review or e.g. 
Ref.Ul). 

An electrolyte solution can be considered as a dilute charge fluid where the closest approach distance between the 
center of a charge with species a and the dielectric wall takes the same value b for all species. The reason is that the 
differences in the various ion diameters are negligible with respect to all other characteristic lengths. (6„ = h for all 
a's, whether 6„ is only determined by the radius of the excluded- volume sphere of species a or 6„ involves some other 
more complicated microscopic mechanism for the short-distance repulsion from the wall. For instance, a layer of water 
molecules, with a thickness of molecular dimensions, may lie between the wall and the electrolyte solution, as it has 

3|.) As a consequence, as 
shown in Ref.[l|, called Paper I in the following, the coefficient /„„' (x, x') of the 1/y^ tail of w„„' (x, x', y) is a product 
of effective dipoles D^{x) and D^i{x'), 



been suggested for instance for another situation, the mercury-aqueous solution interface 






w^^,{x,x\y) ^ ^3^ ^ 2 

(Therefore the tail of Wcc,{x, x' , y) between two particles of the same species a is repulsive when x = x' , as it is the 
case for identical point dipoles with the same direction.) 

The general result ^ arises from a property about the screened potential 4> defined as follows. {SqSq' /eso\v)4> is 
the immersion free energy between two infinitesimal external point charges Sq and Sq' calculated in the framework 
of the linear- response theory as if the radii of the excluded- volume spheres of the fiuid charges were equal to zero pi- 
(The effect of hard cores is briefiy discussed after ll65ll .^ As shown in Paper I, when all particles have the same closest 
approach distance b to the wall, 

y^+co y^ 

(In ^ D^{x) vanishes for x < b.) In the following, a quantity that is independent of charge species a is denoted by 
an overfined letter when it is analogous to another one that depends on a, as it is the case for D^{x) and D^{x). 
Since (f> obeys an "inhomogeneous" Debye equation where the effective screening length depends on the distance x 



from the wall through the density profiles, D^{x) has the same sign at any distance x from the wall, contrarily to 
the effective dipole D^{x), the sign of which may a priori vary with the distance x. Thus, the 1/y^ tail of (j>{x, x',y) 
is repulsive at all distances x and x' from the wall. Moreover, the a;-dependent screening length tends to the Debye 
length £,0 at large distances, and D^{x) can be rewritten, for x > b, as 



2e„ e-''°("=-'') 



cxp , 



D4x) = ~J^^ l + C^ + G^'ix) , (4) 

where G^ (x) tends to zero exponentially fast over a scale of order 1/k-d. In Gauss units, the Debye length ^^ reads 



where p^ is the bulk density of species a. C^ is a constant which vanishes, as well as G^ (x), in the infinite-dilution 
and vanishing-coupling limit considered hereafter. The global minus sign in Q has been introduced, because, in the 
latter limit and in the case of a plain wall {e„ = Csoiv), D^{x) is expected to have the same sign as the dipole d{x) 
carried by the set made of a positive unit charge and its screening cloud repelled from the wall. The sign of 1 + C^ 
depends on the temperature, on the composition of the electrolyte, on the value of the closest approach distance h to 
the wall, and on the dielectric constants e^ and esoiv 

In an electrolyte solution, the Z„'s of all species a's are of unit order and the diameters ajs, of excluded- volume 
spheres also have the same typical value, denoted by a. Moreover, all densities pf s are of the same magnitude 
order. Thus, if the solution is highly diluted, the Coulomb coupling between charges of any species separated by the 
mean interparticle distance a is weak: the condition of low densities, a/a <C 1, implies that (/3e^/(esoiva)) ^ 1, if the 
temperature is high enough for /?e^/(esoivO') to be far smaller than 1 or of unit order. Detailed scaling regimes are given 
in Section ITbI In the corresponding limit, denoted by the superscript (0) hereafter, where the fiuid is infinitely diluted 
and extremely weakly-coupled, the large-distance behavior w^^,(r,r') of the effective pairwise interaction w^^i{r,Y') 
is the same as if the charges e„ and e„/ were infinitesimal external point charges embedded in the infinitely-diluted 
and vanishingly-coupled fluid, 

y;-jO) ^^Z„Z„,0"^(O). (6) 

Moreover, in this limit, the density proflles are uniform at leading order and D^ {x) is given by (@J| where the constant 



C^ and the function Gj, {x) vanish: C\'=Q and Gm {x) 



im 



. Then, by virtue of ^ and ®, 



Di''{x) = Z^D^ {x) with D^{x)^-J . (7) 

Y ^solv V ^solv ^u 



As long as the dilution is high enough, the large-distance behavior w^l^,{r,r') of the effective pairwise interaction 
w„„'(r,r') is expected to have the same functional form as its expression w^^j (r,r') in the infinite-dilution and 
vanishing-coupling limit. This assuption is supported by two reasons. First, since densities are low, the functional 
form of the large-distance behavior w^^,(r,r') is ruled by the effect of long-range Coulomb interactions, whereas 
short-ranged hard-core repulsions are only involved in the values of the coefficients of this leading tail. Second, 
the leading Coulomb effects are due to the large-distance nonintegrability of Coulomb interaction, and the leading- 
order contribution from any integral involving the Boltzman factor of either the effective interaction w„„' or the 
bare Coulomb interaction Ca^e^'V is obtained by Hnearizing the latter exponential factors. This is the procedure that 
underlies the Debye-Hiickel approximation for bulk correlations, which was initially derived as a linearized Poisson- 
Boltzmann theory, where w^^> is dealt with in a linear-response framework as if charges e„ and e^> were infinitesimal 



external charges 



iU 



(In the Mayer diagrammatic approach of Debye-Hiickel theory, the linearized Boltzmann 



factor is that of the bare potential and one must also resum the infinite series of the most divergent integrals that 

n 

arise from this Hnearization (see e.g. Ref.[3|).) The second reason amounts to state that leading Coulomb effects are 
properly described in a linearized mean-field scheme. 

In other words, as long as the dilution and the temperature are high enough, in w^^, the many-body effects beyond 
the linearized mean-field structure only result into the renormalization of charges and of the screening length with 
respect to their values in the infinite-dilution and vanishing-coupHng limit, namely with respect to the bare solvated 
charges ZJs and the Debye screening length ^o- For instance, in the bulk, wf^f(r, r') behaves as the leading-order 
function w^^f (r,r') given by the Debye-Hiickel theory for point charges, but bare charges ZJs and the Debye 
length K~^ are replaced by effective charges Zf °*^'s and the screening length k~^, respectively: 

p2 -KB|r-r'| 

^0ffB^cffB_^ -, (8) 



|r-r'|^+oo esolv 

where 

Zf^(°)=Z„ and 4*'^=Ko. (9) 

(In the following, the superscript B signals all bulk quantities.) In the present paper we show that the effective dipole 
D^ {x) in the large-distance pairwise interaction along a dielectric wall takes the form 

D.{x) = -J J-^ Zf-" [1 + Gl^^{x)] , (10) 

V ^solv Y ^solv /^ 

where G'^^p(x) is an exponentially decaying function which tends to when x goes to infinity. At distances from the 
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wall larger than a few screening lengths, D^{x) takes the same functional form as D^^>{x), and many-body effects 
reduce to the introduction of effective charges Z'^™ and of a screening length k~^. 

On the other hand, in the bulk many-body effects upon effective charges and the screening length arise only from 
pair interactions between ions, whereas along the wall they involve also the electrostatic potential and the geometric 
constraint created by the wall. In order to investigate these differences, we also determine the bulk pairwise interaction 
w^^, up to the same order in the coupling parameter as for the pairwise interaction along the wall. 

B. Main results 

In the present paper we determine wf^/(|r — r'|), D^{x), and D^{x) up to first-order in the dimensionless coupHng 
parameter 

e=-A.,^^«l. (11) 

^ ^solv 

in regimes where e ^ 1, (cr/a)"^ <C 1 and n^^h ^ 1. e oc (/Je^/egoivo) and Ijllll implies that fie^/esow ^ a <C £,ri- 
(Up to the factor 1/2, e coincides with the so-called plasma parameter of an electron gas.) As shown in Section ITTTl 
contributions from steric effects involving a/a are corrections of higher order with respect to the terms of order e in 
some scaling regimes of high dilution where {cr/a)^ <^ e. Moreover, as shown in Section FV C 31 the first corrections 
involving K^b appear only at order ehi{Kob). In the first scaling regime, the density vanishes while the temperature 
is fixed; then /3e^/(esoivO') and l3e^ / {ssoivb) are fixed, namely 

(— j oc e^ and K„b cc s <^ 1 regime (1). (12) 

In the second case, the density vanishes while the temperature goes to infinity, but not too fast in order to ensure that 
{o/aY ^ s\ then both (3 e^ / {eso\v cr) and fie^ / {egowb) vanish. These conditions can be summarized in the following 
way, 

e^ < f-") <£<Kd6<1 regime (2). (13) 

The expressions at leading order in e and K^h in regime (2) can be obtained from those derived in regime (1) for 
a fixed ratio {e/ n-ab) oc Pe^ /{csoivb) by taking the limit where l3e^/{eso\vb) goes to zero while K^b is kept fixed. We 
notice that, when the solvent is water, the Bjerrum length /?e^/esoiv at room temperature is about 7 Angstroem and, 
for concentrations around 10^'' mol/liter, e is of order 10^^ and {cr/a)^ is of order e^ for cr ^ 5 Angstroem. For the 
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sake of conciseness, both regimes (1) and (2) will be referred to as the "weak-coupling " regime and we shall speak 
only in terms of e-expansions. 

Our exact analytical calculations are performed in the framework of resummed Mayer diagrammatics introduced 



in Paper I. For the inverse screening length in the bulk we retrieve 



Kb — l^rt 



, SsX In 3 , , 



that, up to order e, 

(14) 



where 



s^ = £pfz:". (15) 

q;=1 

The leading corrections involved in the notation o{e) are given in ijfiill . As announced above, only screening effects of 
the non-integrable long-range Coulomb interaction are involved up to order e; the diameters (T„'s of charges appear 
only in higher-order terms. This property holds in the bulk as well as in the vicinity of a wall (see e.g. Ref.^). 
More generally, the specific form of the short-distance steric repulsion between charges does not appear in the leading 
correction of order e. The correction of order e in the bulk screening length k~^ vanishes in a charge symmetric 
electrolyte, where species with charge —Z^e has the same density as species with charge Z„e (E3 =0). If the fluid is 
not charge symmetric, the screening length k~^ decreases when the coupling strength increases. 
Our main results are the following. First, we flnd that 



Z^ff^^ Z <^ l + £ 



Z ^1^ + (^^V (^i - — 
"S2 2 1^2/ U 8 



o{e) (16) 



where (1/6) — (1/8) In 3 > 0. In the case of a one-component plasma, formula lfT6ll is reduced to that found in Ref. 
[121 by diagrammatic techniques. (The expression given for a multi-component electrolyte in Ref.|l3l| corresponds 
to another deflnition of the effective charge and does not coincide with our expression lfT6|l .l As in the case of the 
screening length, there is no correction at order e if the composition of the electolyte is charge symmetric. According to 
the diagrammatic origin of this correction, the contribution to Z'f^^ from a screened interaction via one intermediate 
charge has the sign of Z^Yj^ whereas the contribution to Z"^^ from a screened interaction via two intermediate charges 
always increases the effective charge. We notice that the existence of the non-linear term Zf in Z"^^ implies that 
w'^^, cannot be written as Z^ip^i where V'q' would be the total electrostatic potential created at r by the charge Z^ie 
at r' and its screening cloud in the electrolyte. Vc' does not exist beyond the framework of linear-response theory. 

Second, as expected, the screening length in the direction perpendicular to the wall proves to be the same as in 
the bulk, at least up to flrst order in e. Besides, the renormalized charge Z"^^ deflned in IjToll and the renormalized 



charge Z°^^ in the bulk do not coincide. However, up to order e, their ratio is independent of the species a, 

Zf'"^\l+-f^^^+o{e)\zf'' (17) 



with 



''■'^'■i"(S''"<-')'--')-(l 



2 r 



ac(Aci) In 3 



(18) 



(We notice that the notation o(e) in (I17II contains both contributions such as those in Il(i4|l and terms of order e x K^^b .) 
As exhibited by their diagrammatic origins, the various terms in 7^^^ arise both from the nonuniformity of the density 
profiles and from screened interactions via two intermediate charges. These profiles, which have been calculated 

n 

explicitly in the limit of vanishing n^b in Ref.[2;], result from the competition between, on one hand, the screened 
self-energy arising both from the electrostatic response of the wall and its steric deformation of screening clouds, and, 
on the other hand, the profile of the electrostatic potential drop which these two effects induce in the electrolyte. 
More precisely, C[ written in Ijl46|l is the first-order renormalization of the amplitude of D^{x) (see Q), which 
originates from the nonuniformity of the density profiles. (The screened potential (j) appears as an auxihary object in 
the resummed Mayer diagrammatics, and the expressions of C\ and G^ (x, x') are calculated in Section^) C^ 
and its sign depend on the composition of the electrolyte, on the closest approach distance b to the wall (through the 
parameters K^b and /?e^/(esoiv&)), and on the parameter Ad, which charaterizes the difference between the dielectric 
permittivity of the wall and that of the solvent, 

Ael = '"~''^°'\ (19) 

The second term in the r.h.s. of Ijl8|l originates from the renormalization of the screening length and from the 
difference in the contributions from four-body effective interactions in the bulk and along the wall. The three-body 
effective interactions do not contribute to 7*^^^ - so that 7*-^^ is independent of the species a -, because they give the 
same corrections to the amplitudes of w^^, and w^o'. The constant ac(Aci) is written in Hlfi9|l . If esoiv > £»-, as it is 
the case when the solvent is water and the wall is made of glass, ac(Aci) > (1/2) In 3 and the second term decreases 
the ratio Zf'^/Zf''. We notice that CJ'p'^^^x) at first order in e is given in Section lyiDl Contrarily to £»i°)(a;), 
the sign of D^\x) may vary with the distance x, and it depends on the composition of the electrolyte, on the closest 
approach distance to the wall b, and on the ratio between the dielectric constant of the wall and that of the solvent. 



Contents 



The paper is organized as follows. The large-distance behaviors of the effective pairwise interactions w^^i , in the bulk 
or along the wall, are investigated through the large-distance decay of the Ursell function h^^i according to relation 
Q. The latter decay is conveniently studied from Mayer diagrammatics generalized to inhomogeous situations. In 
Section lUl we recall the resummed Mayer diagrammatics introduced in Paper I in order to systematically handle 
with the large-distance nonintegrability of the bare Coulomb potential (far away or near the wall). There appears 
a screened potential cfi, which coincides with the interaction defined from the immersion free energy between two 
infinitesimal external point charges (see Section III B|l . In the bulk, is a solution of the usual Debye equation. Near 
the wall (p obeys an inhomogeneous Debye equation, where the inverse screening length depends on x. In Section Hi ("I 
a decomposition of h^^' into four contributions enables one to show how the basic internal- and external-screening 
sum rules arise in resummed Mayer diagrammatics, and how they are preserved if only some subclass of diagrams is 
retained. It also allows one to show that, for a symmetric electrolyte, J2a P" (i")^°a' (i"i ^') decays faster than h^^i (r, r') 
(see Sections llV El and IVI B|l . We also recall the graphic reorganization of diagrams devised in Paper I for the study 
of the general structure of large-distance tails in dilute regimes. 

Systematic double-expansions in the dimensionless parameters a/a and e can be performed from resummed dia- 
grams. In Section Hill we exhibit the nature of the first various contributions. This leads us to introduce the scaling 
regimes l(T2)l and lfT3JI where the correction proportional to the coupling parameter e is the leading contribution. (We 
also recall the expression of the pair correlation at any distance at leading order in e.) 

Sectionl^is devoted to bulk correlations. We take advantage of the full translational invariance in the bulk in order 
to resum the four geometric series which appear in the Fourier transform of the graphic decomposition of h^^, recalled 
in Section III CI Thus, we obtain a compact formula for the large-distance behavior of /if„/, where the contributions 
of both charges are factorized. This formula is appropriate to obtain systematic e-expansions of the screening length 
and of the renormalized charge from the e-expansions of resummed diagrams. 



In Section IIVDI we also show how to retrieve the corresponding corrections of order e by a more cumbersome 
method which will be useful for the calculations in the vicinity of a wall, where the translational invariance is lost 
in the direction perpendicular to the wall. In position space every convolution in the graphic representation of h^^, 
decays exponentially over the Debye screening length k~^ at large relative distances r, with an amplitude which is 
proportional to 1/r times a polynomial in r. The resummation of the series of the leading tails in r at every order 
in e must be performed in order to get the exponential decay over the screening length k~^ calculated up to order e 
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(see Appendix El- On the contrary, the correction of order e in the renormaUzed charge can be retrieved from only a 
finite number of resummed Mayer diagrams. 

In SectionElwe recall how the screened potential (f>{x , x' , y) and the effective dipole D^{x) in its large-y tail are 
formally expressed in terms of the density profiles in the vicinity of the wall 1|. Then, the e-expansion of D^{x) can 
be performed from the e-expansion of the density profiles, by applying the method devised in Ref.p|. The density 
profiles, which vary rapidly over the Bjerrum length /3e^/esoiv in the vicinity of a dielectric wall, have been expHcitly 
determined up to order e in the limit where Kd& vanishes in Ref.p], and we explicitly calculate D^{x) up to order e 
in the same limit. 

In Section IvTI we recall how the structure of the effective dipole D^{x) in the 1/y^ tail of h^^'{x,x',y) is given in 
terms of the graphic representation written in Section Hi CI We also derive a sum rule for ^^ e„p„(a;)-D„(a;). By 



using the resummation method checked for the bulk situation in Section HV PL we determine the renormalized value 
of the screening length in the direction perpendicular to the wall at first order in e. For that purpose, in Appendix 
El we show that the leading term in x at every order e'' is proportional to {x — b^ exp [—Ko{x — b)], and we resum 
the series of these leading terms. Thus, we check that the correction of order e in the screening length is indeed the 
same in the bulk and in the direction perpendicular to the wall. Then, D^{x) is determined up to order e by only two 
resummed Mayer diagrams. Explicit calculations are performed in the limit where Hub vanishes and the expressions 
of Z°^^ and Z°^^ are compared. Their physical interpretation is given thanks to the diagrammatic origins of the 
various contributions. 

II. GENERAL FORMALISM 
A. Model 

In the primitive model defined above, the hard-core effect between two species a and a' can be taken into account 
in the pair energy by an interaction Vse. which is infinitely repulsive at distances shorter than the sum {a^ + cr^>)/2 
of the sphere radii of both species. Its Boltzmann factor reads 

f if |r-r'| < ((7„+a„0/2, 
exp[-/3z;sfl(|r-r'|;a,a')] = < (20) 

[ 1 if |r-r'| > ((7„+a„0/2. 
Since charges are reduced to points at the centers of excluded- volume spheres with the same dielectric constant as the 
solvent, the Coulomb interaction between two charges can be written in the whole space (even for a; < or a:' < 0) 
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as (Z^Z^ie'^ /eso\v)v{r,r'), where t;(r, r') is the solution of Poisson equation for unit point charges with the adequate 
electrostatic boundary conditions. Since the half-space a; < is occupied b y a material with a dielectric constant ew, 
v{r, r') in Gauss units reads, for a; > and a;' > and for any |r — r'| > 



3 V a n 



Aoi, defined in Ijl9|l . lies between —1 and 1, and r' is the image of the position r' with respect to the plane interface 
between the solution and the dielectric material. In the bulk the Coulomb potential reads 

-«(|r-r'|) = ^. (22) 

The total pair energy C/pair is 

t/pair = -^vsn{\n-rj\;ai,aj) + -^- Z^^Z^^v{ri,rj), (23) 

^ . , . ^ . , . ^solv 

where i is the index of a particle. 

In the vicinity of the wall, one-body potentials appear in the total energy of the system. For every charge a self- 
energy 2'f (e^/esoiv)^cit arises from the work necessary to bring a charge Z^e from x = +oo (in the solvent) to a point 
r in the vicinity of the wall. According to 12111 . the wall electrostatic response is equivalent to the presence of an 
image charge —AdZ^e at point r* inside a wall that would have the same dielectric constant Csoiv as the solvent, and 

l^scif(a;) - -Aci^. (24) 

In the case of a glass wall in contact with water, the relative dielectric constant eviz/csoiv of the wall with respect to the 
solvent is of order (1/80) < 1, Ad defined in iflfljl is negative, and VJscif is a repulsive potential. The impenetrability 
of the wall corresponds to a short-ranged potential Vsr (x) , the Boltzmann factor of which is 

{0 if a; < 6 
(25) 
1 if a; > 6, 

where b is the closest approach distance to the wall for the centers of spherical particles, which is the same for all 
species. The confinement of all particles to the positive- a; region and the electrostatic self-energy may be gathered in 
a one-body potential Vwaii, 

K-aii = Y. ysRix^;a^) + ^ ^Zf ^ F.eit (a;,) . (26) 

^solv 
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B. Generalized resummed Mayer diagrams 

By virtue of definition (QJ, the leading large-distance behavior w'^^, of w„„' is proportional to the large-distance 
behavior h^^^, of /i„„' , 

hll, = -(3wZ>, (27) 

because any power [w^^/]", with n > 2, has a faster decay than w^^,. In an inhomogeneous situation h'^^, is 
conveniently studied by means of the Mayer diagrammatic representation of /i„„' . However, the large-distance behavior 
of the Coulomb pair interaction ^(r, r') is not integrable, and every integral corresponding to a standard Mayer diagram 
that is not sufficiently connected diverges when the volume of the region occupied by the fluid becomes infinite. 

As shown in Paper I, thanks to a generalization of the procedure introduced by Meeron [1^ in order to calculate 
/Iqq' in the bulk, the density-expansion oi h^^i in the vicinity of the wall can be expressed in terms of resummed Mayer 
diagrams with integrable bonds F. Since the procedure for the systematic resummation of Coulomb divergences relies 
on topological considerations, the definitions of Mayer diagrams with resummed bonds are formally the same ones 
in the bulk or near the wall. The two differences between resummed diagrams in the bulk and near the wall are 
the following. First, near the wall the point weights are not constant densities but a:-dependent density profiles. 
Second, the screened potential arising from collective effects described by the systematic resummation of Coulomb 
divergences is no longer the Debye potential, but it obeys an "inhomogeneous" Debye equation. 



Ar(/)(r, r') - K^{x)(t>{r, r') = -47r(5(r - r'). (28) 



In (|2Hl K^(x) is defined as 



K^{x) = 47r/3-^ Y^ ZIp^{x), (29) 

Csolv 

a 

where all densities Pa(x)'s vanish for x < b. (f> obeys the same boundary conditions as the Coulomb potential v: 
(^(r,r') is continuous everywhere and tends to when |r — r'| goes to -l-oo, while its gradient times the dielectric 
constant is continuous at the interface with dielectric walls. We recall that particles are supposed to be made of a 
material with the same dielectric constant as the solvent. 

The two resummed bonds F, called F'^'^ and Fr respectively, are written in terms of the screened potential (j) as 

F^^(n,m) = -^Z„„Z„„,^(r„,r„0 (30) 

^solv 
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and 



Fji{n,m) — exp 



2 



/3e 
-PvsR{\rn - r„,.|) Z„„Z„„Xr„,r,„) 

^solv 



2 



1 + — ^=„^c„0(r„,r,„), (31) 

^solv 



where n and m are point indices in the Mayer diagrams. (In the bond notations, the superscript "cc" stands for "charge- 
charge" and "R" means "resummed". Indeed, F^^ is proportional to the resummed interaction <j>{r,r') between point 
charges; Fn + F'^'^ is equal to the original Mayer bond where the Coulomb pair interaction w(r,r') is replaced by its 
resummed expression ^^(r, r'), while the short-range repulsion is left unchanged.) The resummed Mayer diagrammatics 
of /Iqq' is 

, „ r -/v Us 1 

In Ij32|l the sum runs over all the unlabeled topologically different connected diagrams 11 with two root points (r, a) 
and (r', a') (which are not integrated over) and N internal points (which are integrated over) with A^ = 0, . . . , cxd, and 
which are built according to the following rules. Each pair of points in 11 is linked by at most one bond F, and there 
is no articulation point. (An articulation point is defined by the fact that, if it is taken out of the diagram, then the 
latter is split into two pieces, one of which at least is no longer linked to any root point.) Moreover, in order to avoid 
double counting in the resummation process, diagrams 11 must be built with an "excluded-composition" rule : there 
is no point attached by only two bonds F'^'^ to the rest of the diagram. [JJ F]-^ is the product of the bonds F in the 11 
diagram and Sn is its symmetry factor, i.e., the number of permutations of the internal points r„ that do not change 
this product. Every point has a weight equal to Pa.{x) that is summed over all species. We have used the convention 
that, if N is equal to 0, no /^ drnPa,„ (xn) appears, and (1/S'n) [JI ^]n i^ reduced to F{r, r'). Near the wall, A denotes 
a finite-size region bounded by the wall on the left, whereas, in the bulk, A stands for a finite-size region far away 
from the wall. The screened potential 4> is integrable at large distances. (In the bulk 4> decays exponentially fast in 
all directions; near the wall in the large-distance behavior given in Q), D^{x) has an exponential decay and the l/y^ 
tail is integrable.) As a consequence, 11 diagrams correspond to convergent integrals in the limit where the volume A 
extends to infinity inside the bulk or on the right of the wall. 

C. Graphic reorganization of resummed diagrammatics 

In ft,„„'(ra,r^,) we can distinguish four classes of diagrams by considering whether a single bond F'^^ is attached 
to root point a or to root point a' . (a and a' are short notations for the couple of variables (rQ,a) and {Ya',a'), 
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respectively, which are associated with the root points in a Mayer diagram.) h^^i can be rewritten as the sum 

/i„v = hZ, + h^-, + h2, + h--,, (33) 

where in /i^S both a and a' carry a single bond F'^'^, in h^, {h^^, ) only a (a') is linked to the rest of the diagram by 
a single bond F"^"^, and in ft-~J7 neither a nor a' are linked to the rest of the diagram by only one bond F'^'^. 

1. Screening rules 

A first interest of decomposition l(33jl is that it enables one to derive the basic screening rules (recalled hereafter) 
from the fact that they are already fulfilled by the diagram made of a single bond F'^'^ (because of the corresponding 
sum rules obeyed by the screened potential 0(r, r')). Moreover, since the sum rules are Hnked to the large-distance 
behavior of the charge-charge correlation function, decomposition IJ33II also enables one to show that if some diagrams 
are to be kept for their contributions to h~~,^^ in some dilute regime, then the corresponding diagrams "dressed" with 
F'^'^ bonds in h!]^^, ^ h'^, and h~^i are also to be retained, together with the bond F"^"^, in order to ensure that the 
screening rules are still satisfied. 

The basic screening rules are the following. In a charge fiuid with Coulomb interactions, an internal charge of the 
system, as well as an infinitesimal external charge, are perfectly screened by the fiuid: each charge is surrounded by 
a cloud which carries exactly the opposite charge. These properties can be written in a compact form in terms of the 
charge-charge correlation defined as 

C(r,r') ^ e' J ^Z>„(r)<5(r - r') + J] Z^Z^,p^{r)p^,{r')K^,{r,r') i . (34) 

I a a, a' J 

The internal-screening rule reads 

/'dr^Z„p„(r)/i_,(r,r') = -Z„', (35) 

and, by performing the summation ^^, Z„/p„'(r') x l(35|l . the internal-screening sum rule implies that 

f drC{r,r')^0. (36) 

By virtue of the linear response theory, the external-screening sum rule reads 

^ ■ f dr f dr'i;(ro, r')C(r', r) = 1. (37) 

m 



^solv 



The latter equation, derived for inhomogeneous systems by Carnie and Chan 



y 



is the generalization of the sum- 



rule first settled by Stillinger and Lovett [ig for the second moment of C(r, r') in the homogeneous case (see next 
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p„a..p.,^ A. . consequence . .e i„.e™a, .„ee„i„. ™. ™,e. C3 ,.„■* w.a.eve. .„«i.a„ce ,e.„,a..a«„„ 

may be added to the pure Coulomb interaction w(ro,r') [Q\. 

In the bulk, the translational invariance in all directions implies that sum rules It^ffijl and ^l^7^ are relative respectively 
to the fc = value and to the coefRcient of the fc^-term in the fc-expansion of C^(fc). Both sum rules are summarized 
in the following small-Zc behavior, 

^ ' fe^o 47r/3 ^ ^ 



In the vicinity of a wall, there is translational invariance only in directions parallel to the plane interface, and the 
Carnie and Chan sum rule ll^Tjl takes the form of a dipole sum rule 



irectioi 

r + co r + co p 

I dx I dx' / dyx'C{x,x',y) = — 
Jo Jo J 



esoiv ^3g^ 



47r/3 

As shown in Ref. [l3|, the first moment of C{x, x' , y) is linked to the amplitude /c(a;, x') of the 1/y^ tail of C(x, x' , y), 



dyx'C{x,x',y) = ^^2n dx'fcix,x'). (40) 

£«- Jo 

fc{x,x') coincides with —(3J2aa' e'^Z^Z„>p^{x)p^'{x')fa^'{x,x'), where —(3f^^>{x,x')/y^ is the large-distance behavior 

of haa'{x,x' ,y). Therefore, the moment rule (j39ll can be rewritten as a sum rule for the amplitude /^^'(x,^'), first 

derived in Ref.|2(l|. 

f + CXD /'-t-OO 



/.-I-CX3 n+cxi 

/ dx dx' y e Z^Z^ip^{x)p^i{x')f^^i{x,x') 

Jo Jo 



(41) 



8^2/32 • 

(We notice that there is a misprint in Paper I, where the above sum rule is written in Eq.(4) with a an extra spurious 
coefficient l/egoiv on the r.h.s.) 

Now, we show how the combination of decomposition l(33|l with sum rules obeyed by (j) enables one to derive the 
two basic screening rules. A key ingredient of the derivation is the relations between h"^"^ and h~'^ on one hand, and 
h"^^ and h on the other hand, which arise from their definitions. 

In the bulk, because of the full translational invariance, the latter relations take simple forms in Fourier space. 
They read 

Kl^(k) = Fr^,{k) + ^p;F-^(fc)/^;;„,(fc), (42) 

71 

and 

/^r'(fc) = E<^^'?.(^)^"«'(fc)- (43) 

71 
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(For the sake of clarity, in the present paragraph, we forget superscripts B, except in the densities, in h^^i and in C.) 
On the other hand, by virtue of the explicit expression lf59jl of (J>b, 

Y^ Z^plFZ^k = 0) = -Z^,. (44) 

a 

In other words, the part F'^'^ in h^^, already fulfills the internal-screening sum rule. When relations lll2ll and lli3|l are 
inserted in decomposition ^,VA\ of h^^, , 



hl^,{k) = Fll,{k) + 



Y,p-^F:';^{k)h-i,{k) + h-:,{k) 



Y,PiF:^^Ak)h;-,{k) + h:-{k) 



(45) 



Then property II44|I implies that, in J2a ^"Pq^oc'(^ ~ ^)j the contribution from h"^"^ — F"^"^, given in II42|I . cancels that 
from h^'^ , and the contribution from h'^~ , given in (|i3ll . compensates that from h^^ , so that the internal-screening 
rule is indeed satisfied. 

In the case of the bulk external-screening rule IJ38II . the same mechanism operates when the fc^-term in the small-fc 
expansion of C^{k) is considered. The charge-charge correlation Cf^cc, where h^^, is replaced by i^f^/, fulfills the 
second-moment sum-rule. 

Again, by virtue of Ij44|l . decomposition IJ45II implies that the fc^-term in /i~'^ is canceled by the part of the fc^-term 
in h'^'^ — F'^'^ ~ ^ p^ F^^ {k)h~'^^,{k) that arises from the fc^-term in h~'^ . Similarly, the fc^-term in h~~ is canceled 
by the part of the /c^-term in h'^^ = ^ p^_^F^'^_^{k)h~~^,{k) that arises from the /c^-term in h~~ . Moreover, \A?>\ 
and gll imply that Y.^' ^^'Pl'K^'i^ = 0) = - ^^, Zy p'^M:;_'7 {k = 0), so that the part of the k^-teiui in h"'' - F""" 
that comes from the k^ term in F"^"^ is opposite to the part of the fc^-term in h'^~ that is generated by the fc^-term in 

n 

F"^"^. We notice that the present argument is analogous to that found in Ref.|21J| for an analogous decomposition in a 
quantum charge fiuid. 

In the vicinity of the wall, the derivation of screening rules ll35ll and l(^ also relies on the analog of decomposition 
II45II and on two sum rules derived for (j) in Paper I, namely, if x' > b, 

/ dxK^{x) / dy <f){x, x', y) = in, (47) 



and 



/ dx dx'-R:^{x)K\x')f4x,x')^2^^. (48) 

Jo Jo ^solv 
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The Fourier transform of a function /(y) at wave vector 1 is defined as /(I) = Jdyexp(il • y)/(y). Thanks to 
the translational invariance in the direction y parallel to the plane interface, the relations, which arise from their 
definitions, between h'^'^ and h^"^ on one hand, and h'^" and h on the other hand, take the simple form, 

hll,{x,x',l) = FZ,{x,x',l) (49) 

+ / dxi^p^^{xi)F^';^{x,Xi,l)h~^^,{xi,x',l), 
Jo ^^ 

and 

r + oD 

14711 and Ij48|l imply that F'^'^ saturates the internal sum rule Il35|l and the external sum rule Il41|l . respectively. 

The external-screening sum rule lUTjl in the vicinity of the wall is studied again in Section FVI Bl We show that, 
in the case where all species have the same closest approach distance to the wall, decomposition Ij33|l enables one to 
derive a sum rule fulfilled by the effective dipole amplitude Da{x). 

2. Large- distance tails 

Another interest of decomposition (l33|l is that the large-distance behavior of the Ursell function ft,„„' can be con- 
veniently analyzed from this decomposition, after a suitable reorganization of resummed Mayer diagrams, which has 
been introduced in Paper I. The resummed Mayer diagrammatics ll32)l for /i„„' is reexpressed in terms of "graphs" 
made of two kinds of bonds: the bond F'^'^ and the bond / that is defined as the sum of all subdiagrams that either 
contain no F'^'^ bond or remain connected in a single piece when a bond F'^^ is cut. Fh falls off faster than F^'^ at 
large distances (namely, as [F'^'^]^/2) and the topology of subdiagrams involved in / implies that / decays faster than 
F^^ at large distances in a sufficiently dilute regime. Since the reorganization is purely topological, it is valid for 
correlations in the bulk as well as in the vicinity of the wall. According to the excluded-composition rule obeyed by 
resummed 11 diagrams, the functions in the r.h.s. of Ij33|l are equal to the series represented in FigsQEI respectively, 

hZ,ir,r')=F^^{a,a') + J dr,dr[ ^ p,,(ri)p,;(r;)F-(a, 1)/(1, l')F-(l', a') 

71.71 

+ / dridr[ Y^ p-,^{ri)py^{r[) / drzdra ^ p-,^(r2)p^^(r2) 

71 .71 72,72 

F'='=(a, 1)/(1, l')F"'(l', 2)/(2, 2')F'='^(2', a') + ■ ■ ■ , (51) 
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/i::,(r,r') = j dT,,Y,py{v,,)F--{a,c')I{c' ,a') 

7' 

+ / dYc'^P-,'{Yc') j dridr'i ^ p^^(ri)p^j(r'i) 



7 7l:7i 

i^-(a, 1)/(1, l')F'^\l\ c')I{c', a') + • • • , (52) 



while h '^ is defined in a symmetric way, and 

h:-,{r, r') = I {a, a') + j dvj dv,, ^ p-,{v,)p,,{v,,)I{a, c)F''[c, c')I{c' , a') 

7,7' 

+ / *c / drc'^p-,(rc)py(rc') / rfridr^ ^ p^^(ri)p^;(r;) 

7.7' 7i>7; 

/(a, c)F^^{c, 1)/(1, l')F^'(l', c')/(c', a') + • • ■ • (53) 

In the previous definitions c is a short notation for (re, 7), and i stands for {Yi,ai). 

We notice that, according to previous section, any contribution to / automatically generates a change in /i„„' that 
preserves the two basic screening sum rules. In the bulk, the external-screening rule Ij^i8|l is also retrieved from the 
compact formulae obtained by resummations in the graphic expansion (|5Tll - ll53JI . as shown in Section ITV El 

III. WEAK-COUPLING REGIME 
A. Small parameters 

Now we take into account the fact that in an electrolyte all species have charges and diameters of the same magnitude 
orders e and cr, respectively. Moreover, all bulk densities are comparable, and the typical interparticle distance does 
not depend on species: it is denoted by a. First, we assume that the densities are so low that the volume fraction 
(cr/a)^ of particles is small, 

, 3 



(9 «'■ (^^) 



Our second assumption is that the temperature is high enough for the mean closest approach distance between charges 
of the same sign at temperature T, of order /3e^/esoiv, to be small compared with the mean interparticle distance a. 
In other words, the coulombic coupling parameter T between charges of the fluid is negligible. 



r= -^ ex — <1. (55) 

CsolvO \C,n J 

(The proportionality relation in lf55)l arises from deflnition @.) The high-dilution condition ifKijI implies the weak- 
coupling condition T^ <C 1, if Pe^ / {csoivO') is of order unity or smaller than 1. 
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In fact, conditions H54|l and H55|l can be realized in two different kinds of expansions in the density and temperature 
parameters. In the first situation, the density vanishes at fixed temperature; then the ratio between the pair energy 
at contact and the mean kinetic energy, /3e^/(esoivf), is also fixed, namely 



(-) cxT^ case (1). (56) 



In the second situation the density vanishes while the temperature goes to infinity, so that /3e^/(esoivf) also vanishes, 
namely 



r^ < (-) case (2). (57) 



B. Expansions of resummed diagrams 

The discussion of the T- and (o-/a)-expansions of the integrals associated with resummed IT diagrams is easier if we 
split the bond i^R into two pieces 

^R = ^[^'1'+^HT. (58) 



(The notation Frt refers to the truncation with respect to Fr.) Diagrams built with bonds F'^'^, [F'^'^]'^/2 and F^t 
- and with the same exclusion rule for bonds F'^'^ as in 11 diagrams - will be called 11. The Ursell function /i„„' is 
represented in terms of 11 diagrams by the same formula ^IV2^ as in the case of 11 diagrams. The splitting (I58|l has 
already been used for a classical plasma in the vicinity of a dielectric wall in Ref.J3|, and its use was detailed for 
quantum plasmas in the bulk in Refs.|22Ll2(l|. 

For the sake of simpHcity, the scaHng analysis of diagrams is now discussed in the case of the bulk. The bulk 
screened potential (fiB obeys Ij28|l far away from any boundary, where 'k{x) no longer depends on x and coincides with 
the inverse Debye screening length Kd. Then, ll28ll is reduced to the usual Debye equation, and, since 0b is a function 
of |r — r'l that vanishes when |r — r'| goes to infinity, it is equal to the well-known Debye potential (/)d. 



3 — KD|r— r 



M\r~r'\)^M\r-r'\)^- — . (59) 

|r — r I 

The integrals of the diagrams with a single bond can be calculated explicitly, and their orders in F and a/a are the 
following: 

/"drVfF^'(r,r')-0(r°) (60) 
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and 

y'drVf^[F-]2(r,r')-0(r3/2), (61) 

where 0(r°) and 0(r^'^) denote terms of orders unity and T^'^ , respectively. According to IIH1|I 

I -1-F'='=(r,r')- ^[F'"'f{r,r') if |r-r'| < (o-„+cr„/)/2 
FHT(r,r';a,a')= < (62) 

[ ES 7J. [F'T (r, r') if |r - r'| > (a„ + a^,)/2. 

If we assume, for the sake of simplicity, that all particles have the same diameter a, the expression F^^ of Fj„, in the 

bulk leads to 



/ dr'Fjf^(r,r';a,a')= -^cr^ + ^ir- ^^(7^-271 

J "J ^solv 



V 

27r fPe^Z^Z, " ^ 



"J \ ^solv 



^solv 

[C + ln(3Ko(7)] 



(7 



//3e2Z„Z„A™ (-1)" /pe''Z^Z^,Y , , 

-47r ^- ^^ V , ^ ' -^ + i?r3/2, 63 

where i?p3/2 denotes terms which are of relative order F'^'^ with respect to those written on the r.h.s. of IJ63II . 
Therefore, since p^ is of order 1/a^, at leading order / dr' pf FnT(r, r') is a sum of terms with respective orders 

n\ nV, ^^^ r3, r3ln[p)^l, and r3/f^, (64) 

where /3e^/(esoivO') = r/(CT/a) and the function f{u) = En^(~-'^)"""/[(" + 3)!n] vanishes for u = 0. The last term 
in (IMjl arises from the short-distance behavior of the Boltzmann factor, the explosion of which for oppositely charged 
species is prevented by the cut-off distance a provided by the hard-core repulsion. 

As already noticed in Paper I, the contributions from excluded-volume effects in the primitive model are not involved 
in F'^'^ but they are contained in F^. Indeed, the potential 4> solution of II28II describes resummed interactions between 
point charges at the centers of penetrable spheres, because it corresponds to the integral equation 

0(r, r') = v{y, y')-^ I dv" ^ Zlp^{x")v{v, r")0(r", r'). (65) 

^solv J 

a 

We notice that, in the bulk, for the primitive model again, in a linearized mean-field Poisson-Boltzmann theory 
where excluded-volume spheres are taken into account |J|, an extra Heaviside function 9[\r" — r'| — (a^n -f cr^/)/2)] 
appears in an equation analogous to ll65|l . and the effective interaction between two charges e„ and e„' behaves as 
e^e^i exp{— Ko[?'— (cq +o'c,')/2)]}/ {[1 + '«d(o'q + (To')/2]r} at large relative distances r. The latter interaction is equal 
to e^ea'(j>B{r) up to a steric correction of order (kdc)^ cx T{a/a)^. This is also the case in the so-called DLVO theory 



muM 
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2!tI | for another model where every charge is spread over the surface of the excluded-volume sphere instead of 
being concentrated at its center. In the corresponding effective interaction at large distances, the denominator of the 
steric factor which multiplies exp(— KD'')/r takes the slightly different form [1 + KqIcq + ^qO/^]^- The order T{a/a)'^ 
of this steric correction is one among the contributions listed in (Ifiljl . 

By using the variable change r = v/k-d, it can be shown that, when the number of internal points in a 11 diagram 
increases, then the lowest order in F at which it contributes to various integrals also increases. (See e.g. Refs. 23| or 

y.) 



C. Scaling regimes 

As shown in previous section, in the bulk the leading coupling correction is of order F^'^, and the next correction 
without any steric contribution is of order F^ . The orders of the first corrections induced by steric effects involve a /a 
and F through the combinations written in Ij64|l . 

In the first scaling regime lf56jl . all terms in ijfiiji are of order F^, and the leading correction is indeed provided by 
the correction of order F^' ^ arising only from Coulomb interactions for point charges in the Debye approximation. 



Moreover, we notice that in this regime, where the temperature T is fixed, F'^/^ is proportional to \/J2a Pa^I' the 
density-expansions prove to involve powers of the square root of a linear combination of densities. (The appearance of 
such square roots instead of integer powers in density expansions is an effect of the long range of Coulomb interactions, 
which makes the infinite-dilution and vanishing-coupling limit singular.) 

In the second case l|57ll . /3e^/(esoiva') vanishes, and terms in ijfiill are of orders {<j/a)^ and (c/a)^ times a function 
of /3e^/(esoivO') which tends to zero when /Se'^ / (esowcr) goes to zero. The explicit calculations will be performed in 
a subcase where the leading coupHng correction of order F"^/^ is large compared with all corrections involving steric 
effects. This property is fulfilled if {a/a)^/r^^^ goes to zero, and the corresponding subregime reads 

F^ < f-) <f3/2 subcase (2). (66) 

In place of F, we shall use the parameter e defined in lfTT|l . because the first coupHng correction is of order 

F3/2 oc e. (67) 

(See 1I55I1 and the definition @ of Kq-) In the first scaling regime, relation ll56)l can be written as {(j/a)^ oc e^. Then 
all terms in Ij64|l are of order e^, e^ Ine and e^/(/3e^/(esoivO'), and the whole double expansion in powers of s and a/a 
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proves to be a series in integer powers of e times some possible powers of \ne, at fixed /3e^/(esoivO'). In the second 
regime, condition lf57jl reads e^ <C (o'/a)'^ and the extra condition in (l66|l is (cr/a)"^ ^ e. 

In the following, the so-called "weak-coupling" regime refers to the scaling limits Ijl2|l or lH^^jl . Moreover the term 
"e-expansions" refers to e- and a /a- expansions, as if they were always performed in the scaling regime lfT2)l . 

D. Pair correlation at any distance in the weak-coupling limit 



The scaHng analysis of e-expansions for resummed 11 diagrams (see Section IIIIB|I shows that the £-expansions of 
integrals involving / start at least at relative order e. As a consequence, at any relative distance, the pair correlation 
/i^ j, in the infinite-dilution and vanishing-coupHng Hmit arises only from the sum of 11 diagrams with a single bond 
and where the screened potential (j) is replaced by its leading value (j)'-'^': h]^\ — F'^'^^^> + F^ . (In other words, only 
the graph with one bond F'^'^ in h'^'^ and the graph / in h~~, where / is replaced by Fj^ do contribute at finite distances: 
^ccjo) ^ j^cc(o)^ ;,c-(o) ^ ^-^c^(0) ^ q^ ^^^ ^--(0) ^ ^^^3) [F^^")]' + F^? = 4"'.) At any finite distance |r - r'|. 



/i , reads 



h"-"', (r, r') = e ( |r - r'l - :ii^^i^ ) cxp 



2 



-^Z„Z„,0(o)(r,r') 

^solv 

a(0) 



1. (68) 



In the bulk the inverse screening length k in ll28ll does not depend on x,k = Kd and (J)b ~ <f>B given in lf59|l . where 
0B obeys the Debye equation with the same boundary conditions as the bare Coulomb potential Vb far away from 
any vessel surface. Near the wall, since the density profiles created by interactions depend on the coupHng strength, 
'k^{x) has an e-expansion, and so has (l){x,x\y). In the infinite-dilution and vanishing-coupling limit, k'^{x) tends to 
Kn and 0*^°^ obeys Debye equation with the same boundary conditions as the bare Coulomb potential v, which take 
into account the dielectric response of the wall. 

The large-distance behavior of h^^, at leading order, h^\ , is equal to the large-distance behavior of /i^ j, , namely 

C'°^(r,r') = -^Z„Z„-0(")^^(r,r'). (69) 

^solv 

In other words, since Fr decays only as the square of F'^'^, in the diagrammatic representation h'^J', arises only from 



the diagram with one bond F'^'^, where 4> is replaced by 4>^°' . (The diagram with one bond F'^'^ is called Ha in the 
following and is shown in FigO) Subsequently, the first term in the e-expansion of k is 

k(°^ = Kd. (70) 
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IV. BULK CORRELATIONS 



In the bulk, the Ursell function /i„„' decays exponentially fast in all directions |23|- In the high-dilution and weak- 
coupling regime, the leading tail at large distances is a monotonic exponential decay over the screening length 1/k^ 
(see Ref.J23| for a review), while damping might become oscillatory in regimes with higher densities, as expected from 
various approximate theories (see e.g. Refs.|l3| and J22|). 

The resummed Meeron diagrammatic expansions used in the present paper enable one to retrieve the existence of 
an exponential decay in the dilute regime. Indeed, all resummed Mayer diagrams 11 are built with bonds F'^'^ Ij3()|l and 
Fr ll3T|l , the large-distance decays of which are ruled by the screened potential (pB that is the solution of lf28jl in the 
bulk. By virtue of Ij59|l . (J>b falls off exponentially over the length scale 1/ko defined in ©. The monotonic exponential 
decay of ft-„„' over the length scale 1/kb in the dilute regime is expected to be given by partially resumming the tails 
of n diagrams, which decrease exponentially over the scale 1/kd, though the convergence of the corresponding series 
is not controlled. 

Before going into details, we introduce the following definitions. Let /(r) be a rotationally invariant function that 
decays exponentially fast at large distances r. Let Ko be the smallest inverse screening length in the exponential tails 
of /. / may contain several tails exp[— Kdv] / r'^ with various exponents 7's, which may be negative. We define the 
slowest of the exponential tails of /, denoted by /^^"""{r) hereafter, as the sum of all tails exp{— n^r) / r'' , with any 
exponent 7. In other words, /^^"""{r) is the large-distance behavior with the largest screening length in the exponential 
and all possible powers of r. The notation f^^{r) will be restricted to the leading tail in the large-distance behavior 
of f{r): f^^ir) is the leading term in f^^°^{r), namely the contribution in /*''°^(r) with the smallest exponent 7. 
For instance, as argued in Appendix El if / = 0d * ['/'d]^, /^'°"('') = aexp{—Kor)/r and f^ = ysiow^ whereas, if 
/ = 0D * [(/)d]^ *<Pu, /''°'"(r) = [b + cr]cxp{-Kj,r)/r and /^^ = ccxp(-Kor-). 

A. Resummations of geometric series in Fourier space 

The translational invariance in the bulk impHes that the graph series in the decomposition ll33|l - (l53ll of h„^' are 
sums of convolutions. In Fourier space, they become geometric series which are resummed into compact formulae. 
^QQ^(^) merely reads 

^c..^,^^ ^ ^^Z^Z^ 0o(fc)_ 

esoiv l + Mk)I{k) 
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where 



Mk) 



in 



k^ 



and 



^solv , 

7,7 



h'^J^{k) is reduced to a fraction 



hZf{k) = - 



pe^ 



An 



-= — z„z„,. 



Csoiv k^ + kI+ 47r/(fc) 
The same geometric series appears in the case of ^^~/^ and h^^, ^ with the results 

f3e^ 4n 



hi-r{k) 



Esoiv k"^ + kI+ 47r/(k) 



Z„5]p^,Zv/(fc;7',a') 



and 



K^'''0^)=I{k-a,a') 



/3e2 



An 



fsoiv k'^ + kI+ Anl{k) ^ 
Finally, h^^, (k) takes the half factorized form 

/3e2 An 



7' 



hl^,{k) ^I{k;a,a) 



Esoiv fe2 + k2 + 47r/(A:) 



Z„+^p^Z_,/(A:;7,a) 



7' 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



From this expression we readily get that J2a ^^P^^aa'i^ = 0) = —Z^i, which is another writing of the internal- 
screening sum rule l(35ll in the bulk case. (We also notice that the writing of h^^,{k) in iTtTIi is analogous to Eq.(2.110) 
of Ref.yil, where the authors consider the Feynman diagrammatics for a field theory, with some short-distance 
regularization, which modelizes a charge fluid.) 



B. Large-distance behavior of bulk correlations 

When the bulk Ursell function /i„„/(r,r') is only a function of |r — r'| and decays faster than any inverse power of 
|r — r'l when |r — r'| becomes infinite, its large-distance behavior /i^'^/(r,r') is determined by the general formulas 
recalled in Appendix El As checked at first order in e in next Section, the singular points of I{k;j,-f') in the weak- 
coupling regime are more distant from the real axis in the upper complex half-plane k — k' + ik" than the pole /cq 
of the fraction l/[k'^ + k^ + I{k)] that has the smallest positive imaginary part. Moreover, fco is purely imaginary, 
ko — iK^, and feg is a pole of rank 1. 
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Therefore, in the weak-coupHng regime, the slowest exponential tail of h^^, is a purely exp(— KBr)/r function, as 
well as the slowest exponential tails of h'^J^J , h~^,^ , h'^',^ and h~~, ^ . By inserting the property 



Res 



k^ + kI + Ai:I{k) 



k—kn 



2ko + 47r 



dl{k) 



dk 



H -1 



(78) 



into the general formula l|A2|l applied to / = h^^i given in ifTTIl with fcg — in^, we find that the large-distance behavior 
/i^f'?(r) of h'^^>{r) takes the form 



^L^'(^) 



/3e2 ^- ^- ---^Br- 



^cff B ycff B _ 



^solv 



where in^ is the pole of l/[fc^ + k^ + I{k)] with the positive imaginary part, and 

'l - i{2T:/K^)d7{k)/dkU^^ 



(79) 



(80) 



C. Large-distance tail at order e 



According to the scaling analysis of Section IlII Bl the first term in the e-expansion of I{k) is /'^'(fc) with /^^^ 
j^ccj2^2. / is calculated from the definition (JZU of I{k) with /^^^ in place of /. By using 



^B](k) = yarctanf — 



and 



imH-^'^ 



3 / ^3 






where S„i is defined in Ijl5|l . we get 



4,7'-'(.)..4(|)^^»-n(A), 



f(i) 



f(i), 



/'^■' and / (fc) have a branch point at A: = 2iKr,, while fc^ + k^ + AttI (k) 



27r-(i) s , , s' 
IH — jl [iKoj + o(ej 



(81) 

(82) 

(83) 
has a pole at the value of k equal to 

(84) 



The leading corrections involved in the notation o{e) are given in Ij64|l . The latter pole is closer to the real axis than 
the branch point at fc = i2Ko- Therefore, at first order in s, the singular point in h^^,{k) that is the closest one to 
the real axis in the upper complex half-plane of k is the pole of l/[k'^ + k^ + AttI (k)]. 

The scenario of Section ITV Bl does happen at leading order in e and the large-distance behavior /iff^?(r) up to order 
e takes the form (j79ll where fco = i[Kn + <5kb'] with 



Sk 






(85) 
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According to (I8l-ill we find 

'-^^J^Y'^^ (86) 

We retrieve the formula of R,ef.|l,l| obtained from integral equations. It is reduced to the results obtained by Mitchell 
and Ninham through diagrammatic techniques for the one-component plasma |l3 or for a two-component electrolyte 
2l|. (The formulae for the one-component plasma can be derived from those calculated for a two-component plasma 
with charges e_|_ and e_ and densities p+ and p_ by taking the limit where e_ vanishes while p_ diverges under the 
constraint e_p_ = — e+p+.) The correction Sk^-' vanishes in the case of a 1:1 electrolyte. If the electrolyte is not 
charge symmetric, the expression (|86)l shows that the screening length 1/k^ is a decreasing function of the coupHng 
parameter e at first order in e. 

The bulk effective charge Z°^^ up to order e is calculated by formula ll871| . The explicit result is written in iflfill . 
In view of the discussion of next Section, we rewrite Z^^^ as 

Zf'' = Z4l + Ai^^+oie)]. (87) 

The amplitude ifZfljl of /if ^f (r) up to order e can be rewritten as 

esniv r 11 " " i J 



To our knowledge, the amplitude of h^^?{r) for a multicomponent plasma has not been calculated in the litterature 
previously. In the Hmit of the one-component plasma, it is the same as that found by Mitchell and Ninham in Ref.|l3| 
by a diagrammatic method. 

We notice that the use of ll8?1|l with / in place of / is equivalent to replacing the diagrammatic series of /i^^?, 
h^^,^ and h^^,^ shown in Figs lll3l bv the corresponding series represented in Figs l4l6l As shown in Appendix IbI 
the correction Jkb' to the screening length in h'^'^^, /i'^^^, and h^^^ arises from the whole series in Figs^IEl andEl 
respectively. 

On the contrary, the first corrections to the effective bulk charges can be seen as arising from only a finite number 
of diagrams in Figs^landEl This will be shown in next section. The property rehes on the following rewriting of 

J2 



h^'^fir) = -^ -^^ ( [1 + A(i) -I- A'^V - ^<> r] 6-'^°'^ + o(e)| . (89) 



We point out that l(89ll is vaHd for any distance r. Indeed exp[— (kq + 5n''^^)r\ = (l — (5kb' r) exp(— Kd?*) + 0(e) for 
any distance, whereas exp[— (kd + 5K}^^)r\ = [l — 5k\^^ r + o(e)] exp(— Kd?") only for distances r < L^aayi = ^u/e" with 
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V < 1/2. The condition v < 1/2 ensures that for r < Lmax every n}^ term with n > 2 in the expansion of the 
exponential expl—SK^W] is indeed a correction, of order e"(i^'^) = o(e), with respect to the 0(e^^'^)-term Sn^^r. 

D. Alternative derivation of the bulk large-distance tail at order e 

In view of calculations in the vicinity of a wall, where the infinite series in f ^ h^^^, h'^^^ or h^^^ can no longer 
be resummed in Fourier space, because of the lost of translational invariance in the direction perpendicular to the 
wall, we show how to retrieve the expression (l88|l for h^^ in a less systematic way than the method involving the 
resummed formulae II7H|I and Ij8()|l . (For the sake of simplicity we omit the indices a and a' for charge species in the 
notation /.) 

In the general method of Sections IIV Al and IIVBI we performed Fourier transforms, then we resummed the four 
infinite series / = h"^^^ , h'^^^, h^'^^, or h^^^ that define /if^/ through the graphs shown in Figs lll3l and we got a 
compact formula for h^^,{k) from which we calculated the large-distance behavior h^^f(r) of h^^,{r). Here, on the 
contrary, in each series / we formally calculate the slowest exponential tail /m°™(^) of every graph /,„, with m bonds 
F'^'^ ^ directly in position space by using the residue theorem, and the large-distance behavior f^^ oi f — h'^'^ , h'^~ , 
h^'^, or h^^ is given by the sum (over to) of the slowest tails /m°"(7')'s in each case. (The slowest exponential tail is 
defined in the introduction of Section HVl 'I 

The second procedure is more cumbersome, because the series sums /(fc)'s have a pole of rank 1 at fc = ikb and 
their inverse Fourier transform decay as eiq){—K-Br)/r, whereas each term /m(fc) in the series has a multiple pole of 
rank to, at /c = inu and its inverse Fourier transform behaves as exp(— KDr)/r times a polynomial in r of degree to,— 1. 
We point out that, in the present e-expansion of f^{r) around its exp[— K^rJ/r limit behavior in the infinite-dilution 
and vanishing-coupling limit (where only the bond F'^'^ contributes), for every graph /,„ we must retain the entire 
slowest tail f^""" - namely the entire polynomial in r -, and we only disregard tails exp[—lKDr]/r with I > 2. (See 
example in IJA5II .1 The procedure is legitimate as long as dilution is sufficiently high. Details are given in Appendix 
Eland we give only a summary in the present section. 

1. General structure of the e-expansion of h^^, 

As long as densities are low enough, the graph / decays faster than the bond F"^"^, and, as shown in Appendix IbI 
the slowest tail /m°^('') of fm{r) is equal to exp[— Kor]/r times a polynomial in r of rank m — 1, J2T=o ^m,p'"^- ^s a 
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consequence, the ^-expansion of the large-distance behavior of /if^, reads 



K^fir) = Y. "■" ^p("' "'; ^)' (90) 



r 

p=0 



where the coefHcient Hp{a,a';e) - which is the the sum of the contributions from h'^'^^, h'^^^, h""^^, and h~~^ - 
arises only from the graphs fm with m > p + I in the series representations shown in Figs 1 1131 (We recall that e is 
a short notation for parameters e and a/a in the scaling regimes Ij56|l and II6()|I . as explained at the end of Section 



Moreover, according to the scaling analysis for e-expansions in Section HUB! all coefficients F^.p in the polynomial 
in the slowest tail of every graph fm with m > 1 start at order e"o+™-i^ where tiq = if / = h^^'^, uq — 1 ii f = h^^^ 
or h^'^^ , and no = 2 if / = h^^^. Therefore, the e-expansion of Hp{a,a';e) starts at order e^, and, after reversing 
the summation orders, 



hl^^{r) = 51 ^' II ^^ ^i'^C"' "') (91) 

g=0 p=0 



with 



H('^)(a,a') = n^l,(/ = /»rJ)- (92) 

II91|I displays that, in the e-expansion of ft.fj'f (f) around its exp[— K-^r] / r behavior in the limit where e vanishes, the 
leading tail of h'^^fir) at order e'' behaves as r' times exp[—Kor]/r. Moreover, by virtue of l(92|l . it coincides with 
the first term in the £-expansion of the leading tail in the slowest exponential decay fg^°Y{r) of the graph fq+i with 
(q + l) bonds F"" in / = h^f . 

2. Renormalization of the screening length 

As shown in Appendix^l for each / = h'^'^^ , h"^^^ , h^"^^, or h ^ , the sum f^*{r) of the leading tails at every order 
£"o+g jjj ^ijg £-expansion of f^ir) around its lowest-order e"" exp[—KDr]/r limit can be performed explicitly. (Indeed, 
the coefficient Oq " I'^no+g q£ ^j^g leading r' exp[— Kd?"]/?" tail at order £"«+« in the e-expansion of /^^ coincides with 
the first term in the e-expansion of the coefficient Fq+i^q of the leading r'^ exp[—KDr]/r term in the slowest exponential 
decay /g!fi'('') of the graph fq+i with {q + 1) bonds F"^"^, and the formal expression of -Fq+i.g in terms of / is given in 
IIB7II .1 f^^*{r) proves to be equal to exp[— (5K*r] times f'^'^^""> ^ the value at the first order e"" of the large-distance 
behavior of the graph /i with only one bond F'^'^. As a consequence, the sum h^^*{r) of the leading r"^ exp[—Kj^r]/r 



tails at every order e'' in the e-expansion of h^^f{r) around its lowest-order exp[— Kd'"]/'" limit, namely 



+00 



1r1. 



reads 



/ir?*(r)^^ij(«)(a,a')e' 
9=0 



hi;^?*{r) = F'=''^°\r)e~^''*^ 
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(93) 



(94) 



It arises from the leading e'^r'' e-}q)[— K-^r] / r tails of h'^J^f only. 

Moreover, (5k* coincides with the first-order correction (5kb' to the bulk screening length (see IJHU) calculated from 



the exact procedure of Section IIVB 



(5< = (54''. 



(95) 



Eventually, the resummation of the series of leading tails at every order in e for h'^'^ ^ proves to be a way to retrieve 
the value of (5kb ' • 

3. Diagrams with slowest exponential tails of order e 



As already seen in Section lTlI PL diagram Ila in Figdis the only diagram whose slowest exponential tail (proportional 
to exp(— Ko^)/*^) has an amplitude of order £° x /3e^/esoiv Diagrams whose slowest exponential tails have amplitudes of 
order e are diagrams lib, n^* and He shown in Figs|H|and|ni (which come from the series h^~ , h~^ and h'^^, respectively, 
where / is replaced by /^^\ as shown in Figs^andj^. The contribution of Ilfc to h'^^fir) reads 



^solv 



Z^Z^e- 



EslnSe^'^o'' 



S2 2 



(96) 



while the contribution of He is 

^solv 



z„z„ 



-||)'^--KI)^G-^ 



(97) 



Indeed, the contribution of diagram lib is proportional to the convolution 0d * \<i>i^ calculated in ij A3p . while the 
diagram He involves the convolution (/)d * ['/'d]^ * 0d, whose expression at any distance is given by \k.f\\ . We notice 
that the diagrams Ilf,, Ilh* and He have already been calculated in the case of the electron gas [1^. 

By virtue of ll96|l and l(97ll . the sum of the slowest exponential tails of diagrams Ilf,, n^* and 11^ coincides with the 
expression H89II of h^J^ up to order e, where Sk"^^ is given in Ij86|l and 






(98) 
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with 



-j-(i) _ Eg In 3 



(99) 



and 



^LV^ef'S^ 



1 ln3 
6 ^~8~ 



(100) 



The — (5kb' r term in Ij89|l comes from diagram He. The term Z^Au^ (Z„/ylry ) proportional to Z^ {Z^i) in A)^> comes 
from rib (rife*), whereas the other term, ^Li , arises from He- 

Therefore, the constants A)^> and A]^, are determined by the exponential tails of only three diagrams Ilf,, lib*, 
and li-c- (See the comment after l(89|l for a comparison with the exact method of Section HV CI l Moreover, as a 
consequence of the analysis summarized in Section HVD 21 the coefficient of the r-term in the amplitude of the slowest 
tail of diagram lie with two F^'^ bonds (See FigEjl must coincide with the opposite of the first-order correction (Skb ' 
in the inverse screening length. 

E. Density-density and charge-charge correlation 



By virtue of the bulk local charge neutrality. 



Y.^.Pi^^. 



(101) 



the Fourier transform of the density-density correlation function takes the form 



E'°»/'=''^""'W = E^'^''^(^;"'"') 



Q:,a' 



Csolv fc2 + k2 _|_ 47rJ(fc) 



Y^plY,p^^Z,I(k-^,a) 



(102) 



while the charge-charge structure factor, defined from H34II . 



reads 



C^(fc) = e^ l^plZl + E plplZ^Z^,K^,{k) 



C^(fc) = f^|4 + 4.7(fc) ['^^ + 4.7(fc)]' 



(103) 



(104) 



47r/? 1^'^" ' ""'~'^' P + k2 + 47r/(fc) ^ 

As announced in Subsection III (" ll the expression 110411 of the charge-charge structure factor C^{k) indeed obeys 
the sum rule l(38jl . which summarizes both the internal-screening sum rule l|36)l and the external-screening sum rule 
(13 7p . (If a phase transition gave rise to nonintegrable algebraic tails in I(r) and subsequent non-analytic terms of 
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order fc'' with 77 < in the Fourier transform I{k), then C^{k) would still vanish at fc = 0, but the coefficient of the 
/c^-term would be different from the universal value in ll38ll . as exhibited in the exactly-soluble spherical model of 

Ref-laa-) 

We stress that the fc^-term in the Fourier expansion of C^(fc) is independent of the short-range potential 
fsfldr — r|;a, a'), which must be introduced in three dimensions in order to avoid the collapse under the attraction 

n 

between charges with opposite signs. This property is a consequence of the internal screening rule [6], and it is 
retrieved from the structure of expression Ijl04|l . On the contrary, the fc^-term in X]q a' PaPa'^Za'i^) given in Ijl02|l 
has not any universal value: it depends on the short-distance repulsion in the generic case. However, this is not the 
case for a symmetric 1:1 electrolyte in two-dimensions [2J,IS3) where the pure logarithmic Coulomb interaction needs 
not be regularized at short distances and is scale-invariant. Then, for point charges, scale-invariance arguments lead to 
a value of the dimensionless second-moment of J2a a' PtPZ'^Za'i''^) that depends only on the coupling parameter /3e^. 
We also notice that formulae (110211 and ljl()4|l enable one to retrieve the leading low-density values of the coefficients 
of the k"^ and k'^ terms in the Fourier transforms of the density-density correlation and of the charge-charge structure 
factor derived for a symmetric 1:1 electrolyte in Ref.| 

When there is no charge symmetry in the composition of the electrolyte, the same argument as that used in Section 
IIVBI implies that, according to (110211 and (I104II . the large-distance behaviors of the density-density and charge-charge 
correlations in the high-dilution and weak-coupHng regime are determined by the zero fco = ikb of fc^ + k^ + I{k) : 
they decay over the same screening length as the correlation h^^t. 

In the case of a symmetric electrolyte made of two species with opposite charges +Ze and —Ze and with the same 
radii, J^aPaJ^-, P^^^i-^O^'^^^'^) vanishes by virtue of the local neutrality (jl01|l and of the symmetries (/(k; ++) = 

/(k; ) and /(k; H — ) = /(k; — +)). As a consequence, J2a Pq^cq'('') and ^^ ^, p^p^,/i„„'(r) decay as /(r; a, a') by 

virtue of (I102|l . /(r;a, a') is expected to decay over the length 1/(2kb), by analogy with the infinite-dilution and 
vanishing-coupHng Hmit where it behaves as the diagram [F'^'^]^/2, which falls off over the scale 1/(2ko). Therefore, 
in this peculiar case, the "screening" length of the density-density correlation is expected to be 1/(2kb) at low density, 
in agreement with the result of Ref. 23|- 



,ne a( 

y. 
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V. SCREENED POTENTIAL ALONG THE WALL 

A. Formal expression of the screened potential 

Near the plane dielectric wall located at a; = 0, interactions create density profiles and lf28jl is an inhomogeneous 
Debye equation, where the inverse squared screening length k^ depends on the distance x to the wall. Moreover, 
(f){v,Y') obeys the same boundary conditions as w(r,r') (defined after l|2?l V 

lim --(r, r)= lim --(r, r) (105) 

x^o- Csoiv <JX x^o+ ox 

and 

lim ^(r,r')= lim ^(r,r'), (106) 

x^b- OX x^b+ OX 

since particles are made of a material with the same dielectric constant as the solvent. In order to take advantage of 
the invariance along directions parallel to the wall, we introduce the Fourier transform with respect to the y variable, 
and we write 

cj^ix, x', y) = ^o y 7^ e-'^-^-^y^^{^^x, n^x'- q). (107) 

In the following, the tilde index denotes dimensionless quantities, such as the Fourier transform (I){k-oX , k-^x' ; c\) and 
the dimensionless coordinate x = k^x. 

The solution of the inhomogeneous Debye equation IJ28II requires one to distinguish only three regions: region I for 
a; < 0, region II for < x < 6 and region III for b < x. In regions I and II, k{x) vanishes by virtue of lf29jl . According 
to H28II . the dimensionless Fourier transform (l>{x,x';q) obeys a one-dimensional differential equation. When x' > b it 
reads 

|^-q'W(5?,5?';q)-0 iix<b. (108) 

The solution with boundary conditions Hl()5|l and IJ106II is 

B(iMq|)(l-Aei)el'i|2 ii x<0, 

x,x',\q\)={ ^ (109) 

B{x', |q|) [elil^ - Aoie-lil^] if < x < 5. 



y 



(A similar equation is solved in Ref.^l with a misprint in Eq.(4.20).) 

In region III, when x goes to +00, k{x) tends to the inverse Debye length k~^ ^, and we rewrite the Fourier 
transform of 12811 as 



<^ — ^ - (1 + q^) - U(x) \ (j)(x, x'- q) = -Ati5{x -x') \ix>b 



(110) 
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with 

t/(5?)^i^^Zf[p„(:r)-p-]. (Ill) 

The solution of the one-dimensional equation llllOII can be written in terms of the solutions h of the associated 
"homogeneous" equation (with a zero in place of the Dirac distribution) which is vaHd for — oo < x < +cx3. Indeed, 
the general solution of llllOII for x > b and a;' > 6 is the following sum: a linear combination of two independent 
solutions /i+ and h^ plus a particular solution (/)sing of mi (111, w hich is singular when x — x' and which is calculated 
in terms of h~^ and h~ by the so-called Wronskian method j^. In the following, /i+ (h^) is chosen to be a solution 
which vanishes (diverges) when x tends to +00. In the bulk, k(x) is a constant equal to the inverse Debye length k^: 



/i+ and h can be chosen to be equal to exp[=Fa;-\/l + q^]. When k depends on a;, we look for the solutions h'^ and 
h" in terms of the bulk solutions as 



eT"V^+^[l + H±(J,q2)]. (112) 

Moreover, the particular solutions iJ+ and H~ can be chosen to vanish at x = 6. As shown in Ref. j|, 
and x' > b 



when X > b 



(x, x',q)^ 0si„g(5?, x', q2) + Z(|q|)e-(^+^') Vi+^[1 + H+ {x, q^)] [1 + H+{x', q^)] (113) 



and 



~ 47r I- ~/i / 7 

4i„g(^, x', q) = -— — e-l--- I Vi+q- [1 + H-{mi{x, x'), q^)] [1 + i7+(sup(5?, x'), q^)], (114) 

where inf(a;, x') (sup(x, x')) is the infimum (supremum) of x and x'. Since H^ and H~ vanish at x = &, dH~ /dx\~^-j^ 
is also equal to zero, as can be checked from the formal solutions given in next paragraph. Therefore the Wronskian 
W{q^) takes the simple form 



wic,^)^-2VTT^+'J^:^ 



(115) 

a:— 6 



For the same reasons, the value of ^'dqj) depends only on dH'^/dx\~^i. Indeed, ^(|q|) is entirely determined by the 
ratio of the continuity equations ljl()5|l and ljl()(i|l obeyed by (f) and dcjy/dx at x = 6, and the amplitude B{x' , |q|) in 
region < x < 6 (see 1110911 ) disappears in the latter ratio. 

As shown in Ref.[J|, H^ can be represented by a formal alternating series, which will be used in the following, 

if+(5?,q2) = -r+[l](5?;q2)+r+ [T+[\]] (i;q2)- ... , (116) 
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where the operator T+ acting on a function f{x) reads 



r+[/](J;q2)= / dwe^Vi+^W dte-^^^+^*U{t)f{t). (117) 

Jb Jv 

Similarly H^{x, q^) is equal to the series 

H~{x,<i')^T-[l]{x;q') + T-[T-[l]]{x;<i^) + --- (118) 



with 



T-[f]{x;q^)= rdwe-2\/i+^" [ dt e^^^+^ ' U (t) f (t) . (119) 

Jb Jb 

In the infinite-dilution and vanishing-coupling limit, density profiles become uniform and K^(a;) tends to k^. The 
corresponding screened potential 0^°^ obeys the Debye equation and satisfies the same boundary conditions as the 
bare Coulomb potential v. In other words, H^ and H~ in expressions 111 1311 vanish for 0^°-* (x, x' , q) according to their 



definitions H112|l . while the expression H115|l is reduced to — 2-\/l + q^ for W'^^\q^). Z(°^(|q|) is then determined by 
the continuity equations Hlfl5ll and Hlflfifl . The result reads 

</.(«) (x,x',y) = ^il'.ljr -r') + n^l-^ ,-q.(«oy)^(o)(|q|)g-(2+2')^AW, (120) 

where 



Z(0)(|p|) = 2vr ^^^^;^^ l-A,ie-^'>N(yrT?+|q|)^ 
yrr? (v/TT?+|q|)2-Ae,e-2fc|q| 

The particular solution (Aging (r ~ r') that is singular when r = r' coincides with the bulk screened potential in Debye 
theory, 

0,i„g(r r)- «oy (2^^^^ yiT? 

= 0D(r-r') (122) 

where 0D is written in II59|I . 



B. Large-distance tail of the screened potential 

When X > b and x' > b, 0(x,a;',y) falls off as l/y^, because of the boundary conditions at the interface x — b. 
The reason is the following. The appearance of an l/y^ tail in the large-y behavior of a function /(y) corresponds to 
the existence of a term proportional to |q|, which is not analytical in the Cartesian components of q, in the small-q 
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expansion of /(q) '23|- Functions different from Z{q) in (t)(x,x',q) [see Ijll3|l ] prove to be functions of q^, but the 
boundary conditions at a; = 5 imply that, as well as the small-q expansion of (l){x, x' , q) when x < b (and x' > b) [see 
Hl()9|l ]. the small-q expansion of Z(q) contains a term proportional to |q|. 
As shown in Paper I, the 1/y^ tail of (/) takes the product structure ^ where 



D,p{x) = -\/ ^^ ^[1 + H+{I, q^ = 0)]. (123) 

In H123|l Bz is the coefficient of the |q|-term in the small-q expansion of Z(|q|), 

Z(|q|) = Z{q = 0) + Bz\q\ + 0{\qf). (124) 

We notice that, as shown in Paper I, sum rules obeyed by 4>{x,x' ,y) imply that D^{x) has the same sign for all x's, 
and the l/y'^ tail of <j){x, x' ,y) is repulsive at all distances x and x' from the wall. In H123|l the minus sign in front of 
the square root is a priori arbitrary. It has been introduced, because in the infinite-dilution and vanishing-coupling 
limit and in the case of a plain wall {ew — EsoIv), D^{x) is expected to have the same sign as the dipole d{x) carried 
by the set made of a positive unit charge and its screening cloud repelled from the wall, and H'^{x,o^ = 0) vanishes 
in this limit. 

The large-distance behavior of <f){x^ x' , y) at leading order, cff^ *^°'' , is equal to the leading tail cfr-^' '^ of (j)^^> : cff^ ^^> = 
D^ {x)D^ {x')/y^ with D^ — D^m. D^{a){x) is given by H123|l . where H^ vanishes and Bz is calculated for 
namely, B^ is equal to the coefficient of the |q|-term in the small-q expansion of Z^i^o){q) = Z("'(q). Accordir 
lfT2T| 



Z(")(|q|) = 2^e 



2b 



l-2^^|q| 

^solv 



0(|qn (125) 



and the resulting expression for -0^(0) is written in 10. The expression of the distance yi (x) at which the 1/y^ tail 
in 0^"^ overcomes the exponential tails in (/)(°) has been estimated in Paper I. In the case where the solvent is water 
and where the dielectric wall is made of glass, e»V^soiv ~ 1/80 and y), [x = h) = 7^^, y* [x — b + S,d) — IO^d, 
2/i°^(a; = b + SCz,) = IHo and y^°\x = b + b^o) - 20Co. 



C. Large-distance tail of the screened potential up to order e 

1. Formal e- expansion of the tail 

Because of the nonuniformity of the density profiles in the vicinity of the wall, has an e-expansion. More precisely, 
the £-expansion of the screened potential 4> originating from the e-expansion of density profiles can be determined by 
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111 1311 from the e-expansion of the functions H^ and H~ , which themselves are derived from the formal series Hll(i|l 
and H118|l respectively. 

According to H123II and with the notations of ^ , the first correction D^ {x) in the e-expansion of D^{x) is obtained 
from iJ+(^)(a;, q^ = 0) and from the e-expansion of the coefficient of |q| in the small-q expansion Ijl24|l of Z(|q|), 



B, 



bP +b'z^ +o{e). It reads 



Wi.)^Wi.) 



d!^ + G 



!Xp(l) 



[x] 



(126) 



where the constant Ci is equal to 



C« = 



B 



2B 



(1) 

rjim H+^'\x,e^O), 

^^J a;— i-+oo 



(127) 



exp(l) 



and the function Gj, (x), which vanishes exponentially fast when x goes to infinity, is 



G, 



:xp (1) 



{x) =i?+(i)(i,q2 =0)-^lim H+'^^\x,q^ = 0). 



x—>--\-oc 



(128) 



If we write the e-expansion of ^(|q|) up to order e as ^(|q|) = Z("'(|q|) + Z^^^dqj) + o(e), then B^' = B 



Z(0) 



and i?^ = B^(i). By virtue of H125|l . B^m = — 47r(et,,/esoiv) exp(26). As already mentioned in Section FV Al ^(|q|), 
and subsequently Bz, is entirely determined from the expression of dH^ /dx at x = 6 by the ratio of the continuity 
equations Hlfl5|l and Illfl6|l . When the e-expansions of ^(|q|) and dH^ /dx\i up to order e are introduced in the 
expression 111 1311 of 0(x, x',y), the continuity equations at x = 6 lead to 



Z^'KH) = Z("^(|q|) 



(0)/ 



dH+{x,ci^) 



dx 



(1) 



2yTT^ 



3v/r+?+ |q| - Aeie"2b|ql(3yrTq^ - |q|) 



(129) 



yiT? + |q| - Aeie-2b|ql(yTT^ - |q|) 
Then the coefficient B^w of the |q|-term in the q-expansion of Z(^)(|q|) is determined by using H125|l . and the 
expression of G\ is given by Ijl27|l where 



B 



(1) 



dH+{x,ci^ ^0) 



2B 



(0) 



dx 



(1) 



(130) 



The expression of H~^'^^\koX, q^ = 0) is calculated from '^^ Z^p^{x) as the term of order e in the e-expansion of 
the formal series II116II . The first term in the latter series reads 



T+[l](J;q2)= rdve^-y^^+^ f ^dte-^'^^+^ U {t/n^) _ ^ 

Jb Jv L '^D 



(131) 



As shown in next section, the contribution to k^(x) from each species varies over two length scales, /Se^/cgoiv (times 



Z^) and 1/kd- Therefore, T+[l](a;;q), as well as all other terms in the series II116II . can be expanded in powers of 
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the ratio 2e = Ko/Se^/csoiv of these two lengths. As shown in Ref.jjl, for an operator T+ associated with a function 
similar to k'^{x), the £-expansion of T+[l](5;;q) starts at order e (for any value of b), and the e-expansions of next 
terms in the formal series (IlKill are of larger order in e. Therefore H^^^'> is reduced to the contribution from ^^[1] 

i?+«(J,q2) = -r+[l](2?;q)|('^ (132) 

Similarly, dH'^ /dx\^_~ — — dT^[l\{x; q)/dx\ . We notice that the latter derivative originates both from the deriva- 
tive of T+[l](a;; q)| and of T+[l](x; q)|^ , because the latter term can be written as e times a function of the two 
arguments x and x/s (see AppendixEJ- Similar results hold for H^'-^\ with i?^'^'(a;,q^) — T^[l](a;;q)|'^'. 

2. £-expansion of density profiles 
The density profiles in the vicinity of a dielectric wall have been calculated in the high-dilution and weak-coupling 



regime in Refs. 



y and J3|- (The systematic approach in j| 



is based on the Mayer diagrammatics for the fugacity 



expansions of density profiles. Resummations of Coulomb divergences are performed along a scheme which is similar 
to - but more compHcated than - the procedure used in Section IIIBL because of differences in the topological 
definition of Mayer diagrams in the two cases.) Up to corrections of first order in the coupling parameter s, for Kd& 
and /3e^/(esoiv6) fixed, the density profile reads 






-Zt^^V^''{x;K^) 

^solv 



Pc {x) — p^ exp 

^solv J L \ 

In H133|l 0{e'^) is a short notation for terms of orders written in ifSljl with F oc e'^^^ 



(133) 



More precisely, in Ijl33|l {Z'^e'^ / eso\v)V^^^'^ {x; Kd), called the bulk-screened self-image interaction in the following, is 
the part of the screened self-energy that is reduced to a mere bulk Debye exponential screening of the bare self-image 
interaction (I24II due to the dielectric response of the wall. For two charges separated by a distance 2x, the bulk 
screening factor at leading order is exp(— 2KDa;)- After multipHcation by (3, 



1111 VJi-^y a^i A a •/ liii I ±J- ? 

^solv ^solv ^X y ^solv-^^ 

The other part of the screened self-energy comes from the deformation of the set made by a charge, its screening cloud 
inside the electrolyte, and their images inside the wall, with respect to the spherical symmetry of a charge and its 
screening cloud in the bulk. The deformation stems both from the impenetrability of the wall (steric effect) and from 
the contribution of its electrostatic response if Ad 7^ (polarization effect). When it is multiplied by /3, one gets 

P — Kd L{kdX; Hob) = Zf e L{kdX; K^b) (135) 

Csolv ^ 



38 

^^^' is the electrostatic potential created by the charge-density profile at first order in e. It is given by 



$(1) (x;Ko,Kob,-—-] 
V esoivO/ 



r+oo , r ^e2 



-^ f dx' f dy cj,^'^ (x, x\y)J2 Z^. exp 

^solv Jb J -, 



(1) 



^solv 



{l-Z^^eL{K^x';K^b)} 



(136) 



where (/)^°^(a;, x' , y) is written in (112011 and p-* means that the integral must be calculated at first order in e with Kd& 
and /3e^/(esoiv&) kept fixed. As a consequence, 

Z^ /3e$(i) (x; Ko, K^b, ^^ ) = Z„ e /$ f KdX; Ko&, -^] (137) 

V CsolvO/ V esolvO/ 

L and <i>*^^-' are functions of x which are bounded in the interval < x < +oo, and which decay exponentially fast over 
a few Kq^'s when x goes to +cx). In the case of an electrolyte confined between two walls, the density profile exhibits 
an analogous structure |3?l |. 

3. Explicit results in the limit Kr>b <t^ 1 at fixed /3e /(esoiv&) 

Density profiles have been explicitly calculated at leading order in a double expansion in e and Kj^b with /?e^/(esoiv^) 
fixed in Ref.j2|- Indeed, in regimes where K^b <C 1 the density profile written in H133|l can be explicitly calculated 
at leading order by considering the limit of eL{Kr,x] Kj:,b) and of e/$ (kds; Kd^, /3e^/(esoiv&)) when n^b vanishes at 
fixed /3e^/(esoiv&), and by keeping only the terms of order £ln(Ko&) and e. The corresponding expressions are valid in 
regime (1) where the temperature is fixed (see lll2|l V In regime (2) (see Ijl3|l l. the temperature goes go infinity, and 
the density profiles are obtained from those of regime (1) by taking the limit where /3e^/(esoiv^) vanishes while n^b is 
kept fixed. 

We notice that the corresponding results enables one to calculate the surface tension of the electrolyte- wall interface 
at leading order in e and n^b at fixed I3e^ / (esoivb) c>c e/b I4C|. From the generic expression, one retrieves results already 
known in some special cases. 

In regime (1) Kob vanishes at fixed /3e^/(esoiv&), and the explicit expressions of functions in II133|I are the following. 



L(J; 6) = (1 - A2,) / dt ^ + 0{b) (138) 

A (t + x/F^) -Ael 
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and 



\ esoivoy 






^2 V ^ / ^2 



e-' + ^1^6-2^5^(5?) +0(e5), (139) 



where 0{eb) stands for a term of order eb. By virtue of H13fi|l . the electrostatic potential profile $(^' (x) at first order in 
e arises from the screened self-energy: the term with M comes from the deformation of screening clouds with respect 
to the bulk spherical symmetry, which is described by L H135|l . and the other terms originate from the bulk-screened 
self-image interaction Vi^'"'{x) itTMIl . If A^ = 0, ^(^^(x) is reduced to e {Y.'i/Y.2)M{x). In lfT39|l 



M{x) = / dt- 



2te-^ 



l-A? 



i-(2t)2 ^t^^^r^y 



^el 



C is the Euler constant, 



ff7 = 5 I ^Z, 



^solvb 



where 



giu) = -l 



dt 

t 



S- is defined by 



(140) 



(141) 



(142) 



S±{u 



^3u 



Ei(-3w) ± e"Ei(- 



(143) 



where Ei(— a;) is the Exponential-Integral function: for x > 



Ei(-a;) = - / dt ^C + lnx 



dt 



(144) 



S-{u) decays proportionally to 1/u when u goes to cx), since Ei(— w) behaves as exp{—u)/u for large u. 

We notice that, in the calculation of the part in <i>*^^)(a;) that comes from the bulk- screened image contribution 
Vi^'^^, a key ingredient is the decomposition ljC13|l combined with the expression of the Exponential-Integral function 
114411 ■ g{u) arises because 



dv 



V 



--V9\~ 



(145) 



We point out that g{u = 0) = 0. 

In regime (2) (see ltT3|l l. (ie^ / {esoub) vanishes, whatever the sign of A^ is, because (/Se^/egoiv) ^ & ^ Co- In this 
regime, e vanishes faster than 6, and e/b must be set to zero while 6 is kept fixed in ln6 + gj{e/b); then the latter 



40 



sum is reduced to In b. The result is the same as if the exponential involving the bulk-screened self-image interaction 
V^^^'^{x') in the expression H139|l of ^^^^x) had been linearized at all distances x' , as it is the case in the second 
integral in decomposition ll("13|l . In the following, we write expressions only for the more general regime where e and 
b vanish with their ratio kept fixed. 

In regime (1) (see Ijl2ll l. I3e^ / {csoivb) is finite. For an electrostatically attractive wall (Ad > 0), we cannot consider 
the limit b <C (/3e^/esoiv) ^ ^d, where Aci/3e^/(esoiv^) tends to +oo: there is an irreducible dependence on b. On the 
contrary, for an electrostatically repulsive wall (Ad < 0), we can take the previous limit, where Ad/3e^/(esoiv&) goes 
to — oo. In this limit, b vanishes faster than e, and we must set Kob = at fixed e in the term \nb + gj{e/b) in (|139|l : 
then this term becomes equal to In (|Ad|Z^e/2) + C — 1. 

From the expression of the density profiles up to order e in the vanishing-6 limit, we explicitly calculate the vanishing- 
b limit of the term D^ (x) of order s in the coefficient D^(x). The formal expression of D^ (x) has been derived in 
Section rV(" II The calculations of the zero-6 limits of H^^^>{x^c^ = 0) and dH^ /dx\^ ''{b,c^ = 0) are performed in 
Appendix 10 C\ and G^, (x) in H12fi|l are given by 
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-ln2 + ln& 
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(146a) 




^2 A (i + Vt^-1) -Ad 


■ i+1/2 
2t{t + l) 






(146b) 




i\^2j \Jl (f+Vt2-l)--Ad*' 
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■i + 1/2 
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\ (146c) 
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(146d) 
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(146e) 



and 





Ad S4 f 
' 4 S2I 




S4 r 
^2] I 




1 /S3 

4 VSa 



-{e2SEi(-4x)-Ei(-2i)} 



-2ta: 



dt- 



(i + Vt^^T)- Ad 4^(^ + 1) 



di- 



1-A2, 



1 



Ad S3 
' 6 E2 

Ad /S3 
" 4 VS2 



i + VF^T) -Ad*'-V4 
In 3 



2 



-2ta; 



4t(i+l) 



Y.^P'^Z^g-^ S3 

E^ +E^r ■ 2 



ln6 



1 



"^e2^Ei(-4i) + Ei(-2J) - ^e-"^Ei(-J) - e^Ei(-3x) 



(147a) 
(147b) 

(147c) 

(147d) 

(147e) 



The expressions Ill46a|l and Ijl47all arise from the contribution of the screened self-image interaction V-^'^^'^ Ijl34|l to 
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the density profile ()18;^|l . Thie terms ljl46b|l and ljl47b|l originate from the deformation of screening clouds described 
by the function L given in H138|l . which does not vanish even when e„. = Csoiv The three last ligns in C\ ' /e and 
G^ (^)/£ come from the contribution of the electrostatic potential $(a;) to the density profile. More precisely, 
Ill46c|l - (|147cll . (I146dll - ljl47d|l and H14fiell - H147e|l originate from the functions Af(x), exp(-.x) and exp(-2x)S'_(x), in 
m^fflll . respectively. 

VI. CORRELATIONS ALONG THE WALL 

A. Tails at large distances along the wall 

The Ursell function h^^i cannot decay faster than 1/y'^. Indeed, by an argument based on linear response theory 
and screening in macroscopic electrostatics, the correlation between global surface-charge densities at points separated 



2i|: 



/qq' in the amplitude — /?/„„' of 



by a distance y is shown to decay as l/y^ with a universal negative amplitude 

the 1/y^ tail of h^^i obeys sum rule IjlTll . The latter sum rule holds whether all species have the same closest approach 

distance 6„ to the wall or not. We recall that it is a consequence of external screening, as sum rule (IH8|I . 

On the other hand, as a consequence of the 1/y^ decay of the screened potential cf), according to ij^fijl and l(3T]l . the 
bonds F'^'^ and F^ in resummed Mayer diagrams behave as 1/y'^ and 1/y^, respectively, at large distances y. Since h^^r 
does not fall off faster than 1/y^, no compensation mechanism kills the 1/y^ tail arising from the slowest one among 
the algebraic bonds in the Mayer diagrammatics. Thereore, in a regime where only a finite number of Mayer diagrams 
- or only some infinite class of diagrams - contribute to the large-distance behavior of h^^i , h^^' indeed decays as 
l/y'^. This is the case in the dilute regime studied hereafter. (We notice that if, in some regime, the summation of 
some infinite series of subdiagrams led to an infinite contribution to f^^>{x,x'), then /i„„' would fall off more slowly 
than 1/y^. However, since h^^> is integrable by definition, it cannot decay more slowly than 1/y^.) 

As shown in Paper I, the large-y behavior of /i„„' along the wall is conveniently studied from the decomposition 
described by Ii;^;i|l - lj5^i|l . as in the case of bulk correlations. In the latter graphic representation of h^^', the topology 
of diagrams involved in / implies that the bond / decays algebraically faster than F'^'^ at large distances y (see Section 
III("|I . Moreover, as exhibited in FigsQEl all graphs in h"^^, h~^, h'^~ , and h"" are chain graphs, and, because of the 
translation invariance in the direction parallel to the wall, the chain graphs can be seen as multiple convolutions with 
respect to the variable y. Therefore, every term, except /, in the graphic representation of h'^'^, h^'^ , h'^^ and h 
has l/y^ tails arising from all its F^^ bonds. The 1/y^ tail of every graph in /i'^'^, h'^^ , h^^ and h (see Figs lll3|l is 
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a sum of contributions, each of which is determined by replacing one of the bonds F'^'^ by its 1/y^ behavior at large 
y, while the other part of the graph is replaced by its Fourier transform at the value q = 0. 

Eventually, as shown in Paper I, when all species have the same closest approach distance to the wall. 



h^^'{x,x',y) 



y^+oo 



-P- 



yO 



(148) 



and 



D^[x)^-^{zAd^{x) + C" {x)\+C--[x)\ 

^/esnlv ^ L J J 



Y ^solv 

where C~~(x) and C (x) are related to h^^ and h'^^^,, respectively, by 



(149) 



C:'ix) = fdx"Y,P-,"i^")Z-,"D^ix")h-;;,{x,x",K^q=0) 



(150) 



and 



Z^C" {x) 



I dx"Y,P-,n{x")Z^nD^{x")hl-„{x,x" ,K^(i = G). 



(151) 



An advantage of the resummed Mayer diagrammatic representation is that the contribution from every diagram IT 
can be associated with some physical effect. For instance, diagram Ila made of the single bond F'^'^ describes Coulomb 
screening at leading order, the sum of the two diagrams made of bonds \F'^^^ jl and i^RT respectively contains the 
short-distance repulsion, while diagrams Ilf, lib* and He shown in Figs l8l9l involve many-body corrections to the 
mean-field contribution from 11^. 

In order to trace back the physical effects, we have to identify the contributions to D^{x)D^'{x') from the various 
diagrams 11 defined in Section Fill Bl In other words, we have to recognize in (J148ll - ljl49|l the tails of h'^J^,, /i^~/, h~^, 
and h~~ , the sum of which is equal to h^^' . As shown in the Appendix of Paper I, the latter tails read 



hll,{x,x',y) 






z„z„, 



D^{x) + C" {x) D4x') + C" {x') 



1 



(152) 



^L'(2;'a;''y) 






D4x) + c'^ (x) c'-rix')^, 

yi 



(153) 



h-:.{x,x',y) ^ ^^c:~{x)c:r{x')\. 

y^+oo esoiv y^ 



(154) 
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B. Sum rule for the effective dipole Dc{x) 

According to l|148ll . f^^i{x,x') = D^{x)D^/{x') so that the sum rule l(^ for fcc,'{x,x') can be rewritten as 

j^^ dxY^e^p^{x)D^{x) = ^2%- (155) 

Similarly to what happens for the internal-screening rule l(35|l . the latter external-screening sum rule can be derived 
from the decomposition Ijl49|l . the integral relation Ij5()ll between h'^~ and h~~ , and two sum rules obeyed by 4>, namely 
lUTIl and a sum rule for D^{x) in the f^{x,x')/y^ tail of <j>{x,x',y), 



» /Op 

dxK\x)D4x) = J— ^. (156) 

V £solv 

The latter equation arises from the sum rule ijisjl obeyed by the amplitude /^(a;, x') (derived in Paper I), and from the 
fact that /^(x, x') takes the factorized form D^{x)D^{x') in the case where all species have the same closest approach 
distance h to the wall. 

More precisely, the derivation of Ijl55|l is as follows. The integral relation Ij50|l between h'^~ and h^~ and sum rule 
H47II imply that the contributions from C~^(x) and Z^C (x) to the integral in II155II cancel each other. On the 
other hand, sum rule Ijl56|l ensures that the contribution from Z^D^{x) to the integral in Ijl55|l is already equal to 
the constant in the r.h.s. of the equation. In other words, the bond F^"^ , namely diagram Ha, already fulfulls sum 
rule II155II . 

As a consequence, if some diagrams are to be kept for their contributions to C'^ {x) in some dilute regime, then 
the corresponding diagrams "dressed" with a bond F^^ must also be retained in Z„C {x) in order to ensure that 
screening rule Ijl55|l is still obeyed. 

In the case of a symmetric electrolyte made of two species with opposite charges +Ze and —Ze and with the 
same radii, ^^Pa{x)h^^i{v,v') decays faster than ^-^^'(r, r') in the y-direction, similarly to what happens in the 
bulk (see Section HV E|l . Indeed, symmetries enforce that the local neutrality is satisfied not only in the bulk, 
where p^ = P-, but also in the vicinity of the wall, where p+{x) — p~{x). As a consequence, by virtue of (I149|l . 
p+{x)D+{x) + p^{x)D^{x) = (e/^Csoiv) Yla P<'{^)^'a, (^)- Symmetries also imply that /i^^ = /ill and K^Z = /il+, 
and the definition Ijl5()|l of Cz^{x) yields p+{x)Dj^[x) + p-{x)D^{x) = 0. Subsequently, p+h+^i + p^h^^i decays 
faster than l/y'^. The latter property has been exhibited in Eq.(3.4) of Ref. 2^, where the density-density correlation 
So q' P" {^)p<^' {^')haa' (r, r') in the infinite-dilution and vanishing-coupling limit (where p^ (x) — p^) is shown to decay 
as exp[—2Kjy{x + x')]/y^. 
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£-expansions 



1. Method 



In the general formula l|149|l for D^{x), Z^C (x), as well as C~~{x), is a series of functions, each of which decays 
as a polynomial in x times exp(— Kna^), plus functions which vanish faster (see the general structure Q of D^{x) 
in the introduction). However, since there is no translational invariance in the direction perpendicular to the wall, 
the series C~~{x) and Z^C {x) cannot be expressed as sums of geometric series that could be calculated by such 
a simple formula as ifTTjI . Therefore, the expression of (Z°^^/k) exp[— K(a; — h)] in the large-distance behavior lfT?Hl 
of D^{x) cannot be calculated by the mere determination of the pole of an analytic function and the calculation of a 
residue. 

Though the lost of translational invariance in the direction perpendicular to the wall prevents one from performing 
systematic resummations, D^{x) can be determined up to order e at any distance x (in the sense of the comment 
after Ij89|l ) by the alternative procedure derived for bulk correlations in Section FlV Dl In a first step, the correction 
of order e in the screening length k~^ of the leading exponential decay of D^{x) has to be calculated by the partial 
resummation mechanism whose validity has been checked in the case of bulk correlations (see Section IIV D 2ll . In 
a second step, the amplitude factor in D^{x) up to order e is determined as follows. First, we calculate D^}\x) in 
a form analogous to Ij89|l . which arises from the contributions of only a few diagrams whose amplitude is of order 
e and which decay at large x as exp[— Ko(2: — b)] times a possible linear term in x; in a second step we check that 
the coefficient of the {x — b) exp[— Ko(a; — b)] term, which arises from the second diagram He in h"^"^, indeed coincides 
with the opposite of the first correction to the screening length in the direction perpendicular to the wall, which has 
already been calculated independently. 

The e-expansions of 11 diagrams are more complicated than in the case of the bulk, because the screened potential 
(j) also has an e-expansion when the vicinity of the wall is studied. The first correction to 0'-*'' yields D^ (x) in the 
expression II149II of D^^^x). The leading term in the e-expansion of C~~(x) or Z„C (x) is obtained as follows: 
densities p-,'s are replaced by their bulk values p^'s and both functions (j){x, x',q = 0) and I{x, a;', q = 0) are replaced 
by their leading values (^'"^(x, a;', q = 0) and /(^)(x,.x',q = 0) = {l/2)[F'^'^^°^'^{x,x' ,q ~ 0), respectively. As in the 
bulk case, only the subseries shown in Figs 14161 do contribute to D^^'^{x). 
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2. Renormalization of the screening length 

We recall that, in the bulk case, the leading tail at order e'^ in the e-expansion of the large-distance behavior h'^^J ^ 
of h^Jl^^ around its e-xjp[— K-^r] / r limit decays as r"? times eicp[— n-^r] / r , and the sum h'^fj^^^* of the latter tails decays 
as exp [— (kd + SK%)r], with Sk% ~ (5kb' (see Section lTVD 2^ . In the vicinity of the wall, the contribution from /i^^, to 
Ua (a;) is equal to Z„ Dij,{x) + C (x) according to 1115211 . As shown in AppendixlDl in the e-expansion of Z„C {x), 
at order e"^ the leading term at large x is proportional to {x — b)'^ exp[— Ko{x — b)]. The sum over q of the latter leading 
terms is proportional to exp [—{k-d + 6k*){x — b)] with 

5k*=(5<'. (157) 

According to the discussion of Section HVD 21 the latter partial resummation determines the correction Sk'-^' of order 
e to Ko in the x-direction, Sk* ~ Sk'-^^ According to II157II . the correction Sk'-^' to the screening length in the direction 
perpendicular to the wall coincides at first order in e with the value found for bulk correlations, 

<5«;(i) = <5<>. (158) 

3. Renormalization of the amplitude of D^ (x) 

According to the general method summarized above, the amplitude of D]^' (x) can be determined from only a few 
diagrams in C^ {x) and Z^C (x). Before turning to the expHcit calculations in the regime 6 ^ 1, we interpret 
the various contributions in D^{x)D^r{x')/y^ in terms of diagrams which are representative of physical effects. The 
diagrams that are involved in the determination of the 1/y^ tail of h^^'{x,x',y) up to order in e are the same as in 
the case of the bulk. Diagram Ila in FigO describes the leading screening effect and therefore gives the zeroth-order 
contribution 

Di°)(x)Di"H.') = -^Z^Z^,D^;:\x)D^;:\x'). (159) 

*^solv 

Contrary to the bulk case, because of the nonuniformity of the density profiles in the vicinity of the wall, (j) has an 
e-expansion, and the first correction 0'-^' to 0'-*'' gives a correction of order e in the l/y^ tail of Hq. The contribution 
from diagram Ha to the correction of order e 

[D4x)D^,{x')f^ ^ Df\x)Di]\x') + D'-^\x)D^j\x') (160) 
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reads 



^solv 



-z„z„ 



DTi-)D':\-')+D':\-)D':\-') 



(161) 



The other contributions to 1116011 arise from diagrams lib H;,* and He, shown in Figs l8l9[ where fixed charges interact 
through screened interactions via one or two other charges. (The e-expansions of the contributions from the latter 
diagrams to the l/y'^ tail of h^^^ (x) start at the order e, because they all involve a bond (1/2) [F'^'^] .) The contribution 
to HlfiO|l from diagram !!(, reads 



e y 7t(0), n77— (1). /N 
ZcD^ (x)C„, (X ), 



while lib* leads to a symmetric term in the variables x and x\ and 11^ yields 



(162) 



^solv 



-z„z„ 



c7''\x)d':\x') + d':\x)C7''\x') 



(163) 



D. Explicit results in the limit Kob <S 1 at fixed /3e /(esoivfc) 
1. Separate contributions 

In the limit 6^1, D'/^' {x) is given by Q where the exp(KD&) term disappears in the expression of D^ {x). D'^' (x) 
is calculated from formula H149|l . D^ (x) has been studied in Section and in the limit 5 <C 1 D^ {x) is given by 
II126II . II146|I . and II147|I . The other contributions C^^'-^'>{x) and Z^C"" (x) are obtained by replacing p^"{x") by 
/3^„, (j) by 0^°-', h and h'^^ by the graphs with one / in their series representations, and / by [F'^'^^^>]-^/2 in the 
expressions H15fl|l and H151|l for C^ {x) and Z^C {x), respectively. The contribution from C^ (x) is 



C:-«(x)=<\.x)Z; 



— (1) _oxp(l) ■ 
i>[M + Cxff,! (X) 



>[&] 



'[b] 



(164) 



where 



(1) InSEs 
^[^1 -' — ^, 



(165) 



and 






'[6] 



di- 



^(l-2t)x 



[t + VF~T - Ael(t - A/t^^T)] ' 



2(1 - Ael)2 + 



(l-Aei)^(l-2e-'^)-8Aei(P-l)(l 
t+1/2 



(166) 



with S-{x) defined in 114311 . The contribution from Z^,C (x) reads 



Z^C' ^'\x) ^ Z^Df\x) 



hi3 /S 



x + bI::Ug-^^'\x) 



4 VS2 

In Ijl67|l the coefficient of the Unear term x coincides with — <5kb ' given by IJ86II . while 

1 



B 



(1) 



1 _ ln3 _ ac(Aei) 
3 ~ ^ 2 



S3 

S2 



with 
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(167) 



(168) 



ac(Aoi) 



dt- 



1 [i + vt^^n"- Aei(t-Vi2^n:)]' 



(1 - Ael)2 [I6t4(i + 1/2) - 1] Ael(t - l)(4t2 + 2t + 1) 



8t{t - l/2)2(t + l/2)2(t + 1) t{t - l/2)(t + 1/2)2 



GXp(l) / 



If Csoiv/evv = 80, flc = 1.2, and if esoiv = e^, Oc = 0.84. Gr , '{x) is the exponentially decaying function 



(169) 



^ ^ 4 I S2 



[c] 



dt 



,(1-24)5 



4Ael(t2-l)-(l-Ael)2(i+l) 



[t + Vt^^-Ae,(t-VF-T)]' (t-l/2)(f+ 1/2)2 



+e-^ h5'_(J) + i5'+(S)] + / dt 



-2t. 



2[(l-Ael)2-4Ael(t2-l) 



1 t{t+l){t+l/2)[t + y/W^l-/^,i{t-VW^l)\' 



(170) 



^. Global results 



Eventually, the sum of the various contributions at order e reads 



D^}\x) = 



ln3 /S 



— e- 



("li") 5? + i?(i)+Gr^'n5?) 



^i°H^), 



where B[}'> is the sum of the various constants. 



(171) 



r(i) - r'(i) -X- 7 "r*^^ -X- r(i) 



(172) 



and 



G-p (1) (~) ^ ^exp (1) ^^^ ^ ^^^cxp (1) ^^^ ^ ^exp (1) ^^^ 



(173) 



(^cxp(i)^^-j is a bounded function of order e which decays to zero at least as exp(— a;) when x goes to infinity. The 
coefficient of the term {—x) exp[— Kd^;] in D^y{x) indeed coincides with the first-order correction Sn'-^^ to the screening 
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length in the direction perpendicular to the wall calculated in Appendix ^ with the result II158II . Therefore the e- 
expansion of D^{x) can be rewritten in terms of the explicit expression Q of D^^\x) as 



DM 



2e„ 



X < e 



^solv Y ^solv ^D 

/e-(«B+A-4'*)^ ri + 5a)+Gcxp(i)(j)l +o(£)} 



(174) 



The effective dipole associated with a charge at leading order in e, -D*,°-'(x), is proportional to the mere exponential 
function exp{—K^x). H174|l shows that, when first-order corrections are taken into account, the effective dipole varies 
with the distance from the wall in a more complicated way described by G^^"^ ^^> [x) , the value of which is derived 
from H173|l . The sign of S^^^ + G'™p(^^(2:) may vary with the distance x from the wall and depends drastically upon 
the composition of the electrolyte, the value of the closest approach distance 6, and the relative dielectric constant of 
the wall with respect to that of the solvent. 

Since G^'^p ^^> [x) tends to zero at large x, the effective charge near the wall Z°*^™, defined from the dipolar interaction 
by (Cnj, reads 



r^cffw y 



l + B, 



(1) 



o{e) 



(175) 



where the term 5k^^^ /k^ arises from the \/k coefficient in the definition l|10ll . By virtue of IJ172II and l|86|) 



r^cffw __ y 



l + C-+Z.si^|-e 



-(1) 



ae(Ael) _ 1 

4 6 



|) +o(e) 



(176) 



As exhibited by their diagrammatic origins, the various terms in Z^^^ arise both from the nonuniformity of the 
density profiles described by diagram IIq at order e and from the leading screened interactions via one or two other 
charges that appear in diagrams lib and lib* (Fig|H|l and in diagram He (FigEJ- If Csoiv > £»-, ac( Ad) > (2/3) and 
the four-body effective interactions tend to decrease Z^^™ with respect to its bulk value. 

The comparison of the effective charge Z^^ ™ near the wall with its value Z°* ^ in the bulk given in lIsTIl leads to 

7cffw Ak-'^' 



2-offi 



l + s(l)_^(l) 



+ o(e). 



(177) 



where 



B 



o(i) ^(1). 



W-^(i)=ci^)+i3(V-^ 



(1) 



(178) 



Indeed, according to Ij99|l and Ijl65|l Bry = Aru : the contributions in the bulk and along the wall from diagram Ilf, 
compensate each other and there is no term proportional to Z^ in the ratio Z^^"" /Z^^°. The final result is written 
in eg. 



49 
VII. CONCLUSION 

In the present Paper we have introduced the renormaUzed charge Zf^™ associated with the large-distance dipolar- 
Uke effective interaction between two charges along an insulating wall, when all charges have the same closest approach 
distance to the wall. This charge has been explicitly calculated up to order e in some Hmit of infinite dilution and 
weak coupling when the wall is neutral. 

The renormalized charge could also be calculated in the case of an insulating wall with an external surface charge 
on it, as it is the case for instance when the wall mimicks a cell membrane. Indeed, the general method presently 
devised for the calculation of Z^^"" holds for any density profiles when the e-expansions of the latters are known. On 
the other hand, such £-expansions could be obtained by a generalization of the method presented in Refs.[2| and J3|, 
where the external one-body potential created by the surface charge would be incorporated in the fugacity. 

APPENDIX A: 

The present Appendix is devoted to the determination of the large-distance behaviors of exponentially decaying 
functions, such as those which appear in the resummed Mayer diagrammatics for bulk correlations. When a function 
/(r) is rotationally invariant, its Fourier transform /(k) = /(irexp[— ik • r]/(r) depends only on the modulus k of 
k. On the other hand, when /(r) decays faster than any inverse power law of the modulus r of r, then its Fourier 
transform is an analytic function of the components of k. When both conditions are fulfilled by /(r), the k-expansion 
of /(k) contains only powers of k^. Then the analytic continuation of /(k) = f{k) to negative values of k is an even 
function of k and its inverse Fourier transform can be rewritten as the following integral 

/W=-^-/ dke^'^kfik), (Al) 



where f{k) is a derivable function of k. The one-dimensional integral in ljAl|l can be performed by the method of 
contour integrals in the complex plane k = k' + ik" . (We notice that, when /(r) decays algebraically, then the small-k 
expansion of /(k) contains nonanalytic terms involving either In |k| or odd powers of |k| J37|, and the present method 
does not hold.) 

The slowest exponential tail f^^°^{r) of /, defined at the begining of Section ITVl is determined by the singular point 
of /(fc) that is the closest one to the real axis k" = in the upper complex half-plane. If the latter singular point is 
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a pole fco, its contribution to /(r) is given by the residue theorem 



1 



r'°"M = 2^Res[e^'='-fc/(fc)]|,^,^ 



(A2) 



In the present paper we consider functions f{k) that contain no exponential term and such that ko is purely imaginary, 
ko — iKq. In that case the inverse decay length of f^^°^{r) is equal to the imaginary part of the pole k^. (If there were 
two poles with the same imaginary part and opposite real parts, then f^^°^{r) would be an oscillatory exponential 
tail.) 

Moreover, when /cq = *'^d is a pole of rank 1, /*''°^(r) is a pure exp(— Kd^)/'' function, whereas if fco = *'^d is a 
pole of rank m, /^'°^(r) is equal to exp(— Kor)/r times a polynomial in r of rank m—1. For instance, if / is equal to 
the convolution 0d * [</'d]^, the complete contour integral which determines the r.h.s. of ljAl|l gives the expression of 
<j>-D * [0d]^ at any distance r. According to lf72|l and lIsTjI . (t>r, has a pole at fc = zKd, while [(?!>d]^ has a cut that starts 
at fc = 2iKo and goes along the imaginary axis up to +ioo. As a consequence. 



(r) = 27rln3- 



-Fb{Kor), 



(A3) 



where the first term comes from the residue of exp[ifcr]fc/(fc) at the pole k — iko of (p-oCk), while the second term 
arises from the cut in the definition of [0^](k). At large distances the first term falls off as (po and the second term 
as [(f)'^], because Fb^Kor) decays as a constant times 1/(kd?'). Indeed, 



FbM^-s.J^ '^^ioTTj^ 



-27r [e3''°''Ei(-3Kor) - e'=°''Ei(-Kor)] , 



(A4) 



where Ei(~a;) is the Exponential- Integral function defined in Ijl44|l . (The large-distance behavior of the convolution 
0D * [0d]^ is indeed dominated by the pole of 4>d0<-) at fc = iKo and not by the branch point of the Fourier transform 
of [(?!)d(?')]^ at fc = 2iKD-) In the case of 0d * [4>d]'^, the corresponding slowest exponential tail /'^'°'"(r) is exactly 
proportional to po = exp[~Kor]/r, since the pole at fco = inu is of rank 1. 

For the convolution po * [cpr,]^ * 4>d, [(j)D(k)]^ has a pole of rank 2, and a calculation similar to the previous one gives 



(r) = 167r2 — 



ln3 /ln3 1 

4 ° I 4 3 



where Fc{Kr>r) decays as a constant times 1/(kd''), since 

I (-+00 -2Kurt 

FJK^r) = 327r2— / dt 

>^l Jo [4(1 + t)2 - 



3 — 2kd'' 



KdT 



+ 



■ Fc{Kur), 



(A5) 



(A6) 
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The slowest exponential tail f^^°'^{r) of /(r) — ipo * [4>d]'^ * 0d is given by the pole of f{k) at k — in^, which is of 
order 2, and /^^"""{r) is equal to exp[—Kor]/r times a polynomial in r of rank 1: it is not merely equal to (/)d- 

APPENDIX B: 

In the present Appendix we study the e-expansion of the large-distance behavior h'^l^,. The meaning of e-expansions 
is detailed in Section ITlI CI 

First, we calculate the expression of the slowest exponential tail of the graph fm with exactly m bonds F^'^ in 
the definitions l|5Tll - lf53JI of / — h^^^ , h^^'^, h^^^ or h^^^. (For the sake of simpHcity we omit the indices for 
charge species.) The slowest exponential tail has been defined at the beginning of Section HVl For h'^'^^ and h"^^^, /i 
corresponds to the first graph in Figs Q] and [21 respectively, whereas, in the case of h^^°, /i corresponds to the second 
graph in FiglHl Since F'^'^ decays slower than / (at least in the high-dilution and weak-coupling regime) , 

/i(fc) = k2^.9(fc), (Bl) 

where the poles or branch points of g{k) in the upper complex half-plane are more distant from the real axis than 
the pole k = ikq of l/(k2 + kI). (For h""'', g(k) = 1, while, for ft.'^"^, g{k) is given by ^^ with T(k) = and, for 
h ^ , g{\i) is given by the second term in the r.h.s. of Il76p with /(k) = 0.) Subsequently, the slowest exponential 
tail fi°'"{r) of /i(r) is obtained from IJA2II by calculating the residue of exp[i/cr]/c/i(fc) at /c = zkd, with the result 



/r'°"(o = :^ W'^-khim - ^«o)] L_„ . (B2) 



/f°^(r) takes the form 



-K-D'r 



fiir) = ^—^1.0. (B3) 

According to definitions ll5T | - lf53JI and lf73|l . the Fourier transform of the graph fm with m bonds F'^'^ {rn > 1) reads 

, -, m — 1 



Jm\k) 



-47r/(fc) 



/i(fc). (B4) 



k^ + K^ 

According to the argument leading to l|B2ll . k = iKc is the singular point of fm{k) that is the closest one to the real 
axis in the complex upper-half plane, and the slowest exponential tail f^°^{r) of /m(r) is given by the residue of 
exp[ikr]kfm{k) /{2Tr) at k = inr, according to IIA2|I . Since k = ino is a simple pole for /i, it is a multiple pole of rank 
m for fm, and the latter residue is equal to 



1 1 

X 



27r (m-1)! 



pirn—l 



(B5) 
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The expression ljB5p is equal to exp[— Kq?'] times a polynomial in r of rank m — 1, J^^^o Pin,p'r^, so that 






p=0 



The leading term F^.m-i ^™ '^ exp[— Ko?"] in the residue ljB5|l arises from the (m— 1)'^ derivative of exp[iA:r]. Relation 
ljB4p and comparison of IIB2|I with ljB5|l imply that 



1 



-ZTT 



T m~l 



Fi,o- (B7) 



-^ IIL.IIL — J. / 1 \ I 

(m — Ij! 
The large-distance behavior f^^{r) of /(r) is the sum of the slowest exponential tails of all graphs /„i, 



+ OC 



r{r) = E /™°"(^)- (B8) 



m— 1 

It reads 

+ OC +00 



r{r) = ^-Y.^' E ^™.p' (B9) 

p— m— p+l 



and the large-distance behavior h^^ of h^^, takes the form ijflfijl . 

We now turn to £-expansions of the previous slowest tails. The e-expansion of /(k), written as k^ times a function 
of kJK^ (see l(83ll V generates an e-expansion for /„ through (|Bip and IIB4|I . We recall that /i(fc) is the generic 
notation for the Fourier transform of the graphs with only one bond F'^'^ in Figs 1 1131 Since the e-expansion of I{k) 
starts at order e (see ll83ll V the e-expansion of fi{k) begins at order e"" with no = if / = /i'^'^^, no = 1 if / = h'^^^ 
or / = /i^'^^, and no = 2 if / = h^^^. For the same reason, the £-expansion of fm starts at the order ^"o+m-i^ ^j^j 
so does the e-expansion of the coefficient F^.p of r^ in IIBfill . 

-t-oo 

F™,P = e"° E ^iT'^s'- (BIO) 

q— 7n — 1 

As a consequence, the large-distance tail f^^{r) of /(r) defined in IIB8II has an e-expansion of the form 

/^''(r) = e"°E^''E"p"°^^'' 

p— g— p 

+ 00 q 

.e"o^£<z^a(^«o+9)^P^ (Bll) 



-kdI" 



r 

g=0 p=0 



with ap " '^' = J2m=p+i ^rn°P "^ ' ^nd the e-expansion of hl^{r) has the structure (jHU- In ljBll|l the coefficient 
Qq "^ of £"o+?r'? has a simple expression, 

a(po + q) ^ pino+q) ^ (B12) 
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The leading tail at order £""+« in the e-expansion of f^ is proportional to r' exp[— Kq?']/?'- It comes from the leading 
r'^-term in the slowest tail /^^"(r) of the graph with g + 1 bonds calculated at its lowest order in e, namely e"""*"^. 
The sum f^ * (r) of the leading tails at every order in e in the e-expansion of f^^ reads 

f^*{r) = e"" Y^ «("«+«) {erf . (B13) 

It can be calculated expHcitly by virtue of IIB12|I . because F^-^^^ is given by ljB7|l . where I{k) is replaced by / (fc), 
while Fifl is replaced by the first term in its e-expansion, namely F\^q . We get 






_2^/Vo) 



n 9 



e'^'>F['^°\ (B14) 



|IB14|I implies that the coefficient of r'> in the definition IIB13|l of f^*{r) is indeed such that f'^^*{r) coincides with 
the series of an exponential, 

/-*(r) = -e-('^°+*"B)'-£»oj^^("^o) ^.^jj fc* ^ 2tt ^'"""^ (B15) 

According to l(85ll . Sk^ coincides with the correction of order e in k^, Sk^ — Sk^K By comparison with IIB3|l . the 
relation ljB15p can be rewritten as 

/^''*(r) = e-^"B'- /f ("°)(r), (B16) 

where f^^^^^°' = jj5°*("o) ^^ proportional to ex.p[— Kdv] / r . In other words, the sum f^^*{r) of the leading tails at 
every order in e in the e-expansion of f^{r) around its ex.p[—Kor]/r behavior in the vanishing-e limit is equal to 
exp[— (5kb V] times the large-distance behavior of the graph /i with only one bond F^^ calculated at the first order 
e"». When / = h""'', no = and Z^^'^"") coincides with the diagram Ua = F'"' shown in Figd When / = h"-"" 
or h^^^ , TiQ — \ and f^^^°' is the exp[—Kor]/r tail of diagram Ilf, or lib*, respectively (see FigJHl, the amplitude of 
which is of order e with respect to that of F'^'^. When / = h^^^, uq = 2 and f^^"^'> is the eiq)[~Kor]/r tail of the 
diagram built with (1/2) [F'^'^] * F'^'^ * (1/2) [F'^'^] , the amplitude of which is of order e^ with respect to that of F'^'^. 

APPENDIX C: 

In the present Appendix, we consider the Hmit 6 ^ 1 at fixed Pe^/{esoivb) and we calculate the expHcit values of 
C\ and G^ (x) [see definition 1112611 ] up to terms of order e times In 6 plus a function of /3e^/(esoiv&)- According 
to H127|l , II128II and H13r)|l , the values are determined from i?+ (a;, q-^ = 0) and from its derivative with respect to x at 
point X = b. 
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First, we calculate H~^^^\x,q^ — 0). According to Hi;^2ll . at first order in e, H+(^) = — T+[l]| . In the following, 
the definition II131II of T+[l] is rewritten for q^ = 0, thanks to an integration by parts, as 



T+ [1] (x; q^ = 0) = K{b) - K{x) 



(CI) 



with 



K(x) = 



du 



= -2(«-5) 



Y.pi^l 



Pi 



(C2) 



With these definitions 



X — *- + oo 



(C3) 



and II128II is rewritten as 



-^xp(l) 



(^)=K(i)(J). 



(C4) 



From now on, we consider the regime & ^ 1 at fixed /9e^/(esoiv&)- K^^\x) is determined from the density profiles 
up to order e given in (^Hl- At leading order in the limit where b vanishes, by virtue of ljl^-{8|l and (li;ffl|l . 



^ZleL{x; b) ~ Z^[3e^^^^ { x; n^.b, 






= ei?„(x)+0(e&). 



(C5) 



where Ra{x) is a linear combination of functions of x, where one coefficient involves In 6 plus a function of f3e^/{esoivb) 
in such a way that the limit of this sum is finite if Ad < and 5 = 0. 0{eb) is a short notation for terms of order eb. 
Then, the expression of the density profiles at order e can be rewritten at leading order in b as 



p„(x) 



exp 



^solv 



1 > +eR^{x) + \ exp 



^solv 



■Mm^^(^) 



l}eR^{x) + Oieb,e'^), (C6) 



where 0{e^) stands both for terms of orders written in II64|I with F ex e^'^ (as in Ijl33ll l. and for terms of order e^ 
times a possible sum of a ln5 term and a function of (/3e^/esoiv), which is similar to the coefficient in R^ix) (see 
the comment after (IC5|l V eR^{x) is a bounded integrable function of only x — KqX, while {Pe'^/eso\v)V^^^'^{x) 
is a function of both x/(/3e^/esoiv) and k^x. As already noticed in Section 3.3 of Ref.J2|, an integral where 
{exp [— Zf (/3e^/esoiv)^m'"^(2;)] — 1} is multiplied by e times a bounded integrable function of x is of order e^ times a 
function of 6 and /3e^/(esoiv5), which has a structure similar to Ra{x) in (l("6|l . Thus, according to lj(J2ll and IjC^fill . the 
expression of K{x) at order e, K'^^^x), can be written at leading order in b as the sum of only two contributions 



if (i)(5f) = Kl^Jix) + eKnix) + 0{eb), 



(C7) 



55 
where Ki^{x) and Kr{x) are defined as K{x) by replacing [Pc.{u/ n-o)/ p^ — 1 in (|("2|l by 
{exp [— Zf (/3e^/esoiv)V^m'''^(w/KD)] — l} and i?a(u), respectively. Rc{u) defined in (IC5|I is given by the explicit 
expressions Ijl^i8|l and (ll^fflH . K^ix) is calculated by reversing the order of the integrations ^ du from the definition 
of Kji and Jdi from the expressions of L and $^^^ in R^{u/kj^). Eventually, 



Kb{x) = Kj^{x) + Kj^ix) + KE^pix) + Ks_ (x), 



(C8) 



where the four contributions arising from the terms in the density profiles involving either L{x), M{x), exp(— 5;) or 
S- are written in ljl47b|l . Ijl47c|l . Ijl47d|l and Ill47ell . respectively, by virtue of lj("4ll . 
Now we show that 



Kl'Jix) = Kli^ix)+Oieb), 



(C9) 



where Kll^{x) is deduced from Kiy^ix), defined after IIC7ll . by linearizing the exponential that contains the bulk- 
screened self-image interaction Zf (e^/esoiv)V^,n^'^- Indeed, Ki„-^{x) — Kll^{x) = Q{x) with 



1 /■+°° 
Q{x) = ^ / du 



l-e-2(«-2)^pfZ, 



exp , ^Z^'-e-A 1 - ^Z^'-e~'- 
^ ' 2 °-u / 2 °u 



(CIO) 



Since the functions in the square brackets are positive and 1 — exp[— 2(w — x)] <1 — exp[— 2(u — h)] for x > b, 

^solvb 



0<Q(x)<Qib)^Q[b,s; 



(Cll) 



The double expansion of (5(6, e;/3e^/(esoiv^y) in powers of e and 6 at fixed (ie^ / {tsoivb) can be calculated thanks to 
the following formula (already used in Ref.[3|). We set £„ = Z'^e/2. The function / in the integrand of (|("10ll is a 
function of u that depends on e„ as if / were a function of the two independent variables u and ui — u/e^. We write 
it as 



/ M, — =5(e„ui,ui). 



(C12) 



Since 6 <C 1 the integral jj can be split into the sum of integrals j^ and J~ with b < I and £„ ^ Z <C 1. Then 



Exp 



+ OC 

dufiu,- 
b V £o 



(C13) 



Exp So, I dui Exp g{eaUi,ui) + Exp / du Exp f [u, — 



(//ea)-^+oo 

where Exp denotes an 5-expansion. The identity holds, because when ui < l/e^ then e^ui <C 1, and when u > I 
then (Ec/w) <C 1. When ljC13ll is applied to the calculation of Q{b, e; Pe^ / (esowb)) , the second integral in ljC13|l gives 
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a term of order e^ oc e^, while Exp < 1 — exp — 2(e„Mi ~ b > behaves as —26, so that the first integral provides a 

e„^0 '- L J J 



contribution which starts at order e6; more precisely, 



«i'.^^S)-^^^i?^-'^='' 



(C14) 



where ga has been defined in H141|I - (I145|I . and ©(e^) is equal to e^ times a function with a structure similar to that 
of i?„ (a;) in (l('5|l . The result ()('14|l combined with inequalities lj("ll|l leads to ll('9|l . Eventually ll('7|l can be written 



as 



X(1)(S?) = Kl^ix) + eKnix) + Oieb). 



(C15) 



Kb.{x) is given in (IC8II and, according to its definition. 



K^{^) ^^^^ [e2SEi(-45?) - Ei(-22?)' 

4 2j2 L 



(C16) 



where Ei(M) is the Exponential- Integral function defined in H144|l . 

The derivative dH'^{x,q^ — 0)/dx at first order in s must be performed more carefully. The reason that leads to 

H+^^'l = -T+(i)[l] also implies that 



dH+ 



(1) 



and, similarly to l|C7ll . 





dK 


dx 


The decomnosition l|C13ll lea.ds to 


dKim 


(1) 


dx 


2=b 



(1) 



dx 



dK. 



dK 
dx 



(1) 



dx 



(1) 



+ e 



dKR 

dx 



0{eb). 



2 '^ I Sa S2 



C + ln 



M 



(C17) 



(C18) 



(C19) 



We notice that 



dKu 



dx 






dx 



(C20) 



The reason is that, though Kl^ is only a function of x, K^^ involves a contribution that is equal to e^ times a function 
of the two variables x and x/e, and the derivative of the latter contribution with respect to the second argument x/s is 
of order e. The existence of such a contribution in K^^ is due to the fact that the function Z'^{Pe'^ / esoiv)V^^-^^'^ {x; Kb) in 
the expression H134|l of the density profile varies both over the Bjerrum length /?e^/esoiv and over the screening length 
^n- (This structure arises directly when the equation obeyed by H^ is solved by a multi-scale expansion method.) 
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APPENDIX D: 

In the present Appendix we consider the large-x behavior D'^'^{x) of the dipole D'^{x) that appears in the large-y 
tail D'^{x)D'^,{x')/y^ of h'^Jl^. Calculations are not as straightforwrd as for h'^J^ in the bulk, and we calculate only the 
sum Dl'^*{x) of the leading large-x terms at every order e'' in the e-expansion of £'^^''(a;) around its infinite-dilution 
and vanishing-coupling limit D'^^''^>{x). According to II152|I . 

Dl{x) = -^=Z^ \d^{x) + C''{A . (Dl) 

After insertion of the graphic representation H52|l of If^ in the definition II151|1 of C (x), the latter can be written 

as 

+00 

771—1 

where Jm{x) is the contribution to C (x) from the graph in /i^~„ with m bonds F'^'^(see Figj^J. 

The graph with m bonds F'^'^ in h'^'^,, also contains m bonds /. Therefore, according to the scaling analysis of 
Section lillBL the e-expansion of Jm{x) starts at order e™, Jm{x) = X]q^ Jm {x), where Jm (x) denotes the term of 
order e"^. Therefore, the leading tail at order e* in the large- a; behavior C (x) of C (x) can arise only from the 
leading tails of the Jm(x)'s with m < q. Though we are not able to derive the leading tail of Jm{x) systematically, 
we expect, by analogy with the F^'s in the bulk case, that the leading tail of J,n{x) has the same x-dependence as 
the leading tail of the first term Jm [x) in the e-expansion of Jm{x). As shown herefater, the leading tail Jm ^^{x) 
of Jm (x) is proportional to c'"(x — 6)™exp[— (a; — b)]. As a consequence, the leading tail at order e* in the large- a; 
behavior C {x) of C (x) coincides with the leading tail Jq ^{x) of Jq (x), and the sum C (x) of the leading 

tails at every order s'' in C (x) is C (x) — J2q^ Jq "{x). Similarly, for D1{x) defined in IIDlll 

+00 



Y^solv 



9=1 



(D3) 



The term Jq'{x) of order s'' in the e-expansion of Jq{x) is obtained by replacing every bond F'^'^ by its zeroth- 
order expresssion F'^'^t"), every / by its lowest-order value [F'^'^(°)]^/2, and every weight pdxn) by its bulk value p'^. 
Inspection of Jq (x) for small values of q shows that only the part ^^J = 0d of cfi^"^ does contribute to the leading 
tail Jq ^^{x). Let us denote by Jq{x) the corresponding part in Jq (x). For the sake of simplicity, we relabel 
point pairs {p,p'} = {(rp,7j,), (rp,7p)}, with p = 1, . . . ,q, Xq = Xc' and x'^ = x" , in the opposite sense and we set 
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Up = Xq_(^p_i) and u' = x' ,_^y By using 0, (jl22|l and the change of variable t — 2^/l + q^, we get 

*''^'-^/s (-4(1)7 r '-t r ''"^-"-'••' r '<'-''-" '~' 

2 ^1 Jfc Jb 

Next steps of the calculations involve the following formulae, 

r + oo _ ny. _ 1 



6 ti - 1 ii - 1 



and 



+ 00 



9f -^ ~ -^ 

dui e-l"^-"^l/(?xi) = ^^(4 - 6)e-("^-^) + i?o("2 " &), (D6) 

tj — 1 

where Rp{u) denotes a function whose slowest exponential tail is equal to exp(— m) times a polynomial of rank p in 
the variable x. More generally, we find 

/ d<e-*-l"''-".l(u;-6)''-ie-K-^) = ^l^(^p_&)P-ie-("--'')+i?p_2K-&) (D7) 

Jb tp— '^ 

and 



Jb P 

Eventually, the multiple integral in IJD4II is equal to 

where J^"^ dt2/(t^ - 1) = In 3, and 



(D9) 



Therefore, I?^^''*(x), given by ljD3|l with D^ (x) written in Q, proves to be the series of an exponential whose 
argument is proportional to e{x — b), 

Dl^*ix) = —^ Zj=^- (Dll) 

Y^solv V ^solv /^D 



with 
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By virtue of H8(i|l . Sk* coincides with the first-order correction Sug' to the screening length in the bulk. We notice 
that IDllll can be rewritten as 



Dl'''*{x) = D^°\x) e-'5«*(^-^). (D13) 

IJD13II corresponds to the relation l(Mjl in the bulk case. 
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FIGURE CAPTIONS 

Fig.l Representation of /i^'j^,(r, r') as the graph series defined in Ij51|l . A wavy fine represents a bond F'^'^ and a grey 
disk stands for a bond /. A couple of variables (ri, 7^) is associated with every circle. For a white circle a = (r, a) [or 
a' — (r'j a')], r and a are fixed, whereas, for a black circle i = (r^, 7^), r^ and 7^ are integrated over with the measure 

Fig. 2 Graphic representation of definition II52|I for h'^'^,{r,r'). 
Fig. 3 Graphic representation of definition l|53|l for h^^,{r,r'). 

Fig. 4 Diagrams in h'^^,(r,r') that contribute to the correction of order e in the screening length. A double wavy 
line denotes a bond (1/2) [F'^'f . 

Fig. 5 Diagrams in /i^^, (r, r') that contribute to the correction of order s in the screening length. 

Fig. 6 Diagrams in h^^,{r, r') that contribute to the correction of order e in the screening length. 

Fig. 7 Diagram Ha- 

Fig. 8 Diagrams lib (on the left) and Hb* (on the right). 

Fig. 9 Diagram He- 
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